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Abstract 

In this paper wc consider the circular restricted three body problem which models the 
motion of a masless body under the influence of the Newtionan gravitational force caused 
by two other bodies, the primaries, which move along cicular planar Keplerian orbits. In a 
suitable system of coordinates, this system has two degrees of freedom and the conserved 
energy is usually called the Jacobi constant. In 1980, J. Llibre and C. Simo |LS80b| proved 
the existence of oscillatory motions for the restricted planar circular three body problem, 
that is, of orbits which leave every bounded region but which return infinitely often to some 
fixed bounded region. To prove their existence they had to assume the ratio between the 
masses of the two primaries to be exponentially small with respect to the Jacobi constant. 
In the present work, we generalize their work proving the existence of oscillatory motions for 
any value of the mass ratio. 

To obtain such motions, we show that, for any value of the mass ratio and for big values 
of the Jacobi constant, there exist transversal intersections between the stable and unstable 
manifolds of infinity which guarantee the existence of a symbolic dynamics that creates the 
oscillatory orbits. The main achievement is to rigorously prove the existence of these orbits 
without assuming the mass ratio small since then this transversality can not be checked by 
means of classical perturbation theory respect to the mass ratio. Since our method is valid 
for all values of mass ratio, we are able to detect a curve in the parameter space, formed 
by the mass ratio and the Jacobi constant, where cubic homoclinic tangencies between the 
invariant manifolds of infinity appear. 
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1 Introduction 

The restricted three body problem models the motion of three bodies, one of them massless, 
under the Newtonian gravitational force. Since one of the bodies is massless, it does not cause 
any influence on the other two, the primaries. Thus, the two primaries are governed by the 
classical Keplerian two body problem. Let us assume that these two bodies, whose motion lies 
on a plane, perform circular orbits and that the massless body moves in the same plane. This 
problem is known as the restricted planar circular three body problem (RPC3BP from now 
on). If one normalizes the total mass of the system to be one, the RPC3BP depends on one 
parameter fi which measures the quotient between the masses of the two primaries and that 
therefore satisfies E [0, 1/2]. Then, taking the appropriate units, the RPC3BP is Hamiltonian 
with respect to 

H[q,p, t; /i) = — II — ■ --|| - II . .,.|| (1) 

2 k + MoW \\q - [1 - lJ.)qo{t)\\ 
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where q,p G and —fiqo{t) and (1 — fj,)qo{t), qo{t) = (cos i, sin i), are the position of the 
primaries. This Hamiltonian has two and a half degrees of freedom. It has a conserved quantity, 
usually called Jacobi constant, defined as 



J{q,P,t^,^J') = H{q,p,t;fi) - {qiP2 - q2Pi)- 



(2) 



The purpose of this paper is to study the existence of oscillatory orbits for the RPC3BP, 
that is, orbits such that 



When /i = 0, the motion of the massless body is only influenced by one of the primaries and 
therefore it satisfies Kepler's laws. In particular, oscillatory motions cannot exist. In this paper 
we show that oscillatory orbits do exist for any value of mass ratio ji S (0, 1/2]. 

The study of oscillatory motions was started by Chazy (see |AKN88j ). In 1922, he gave a 
complete classification of all possible states that a three body problem can approach as time 
tends to infinity (see Section 2.4 of |AKN88j ). For the restricted three body problem (either 
planar or spatial, circular or elliptic) the possible final states are reduced to four: 

• (hyperbolic): ||(7(i)|| —5- oo and ||(7(t)|| — c > as t — ±00. 

• (parabolic): ||(7(i)|| —5- 00 and ||^(t)|| — )■ as t — >• ±cx). 

• (bounded): lim supj^j.^^ \\q\\ < +00. 

• OS^ (oscillatory): limsupj_^^oo llg'H = +00 and liminfj_!.±oo Ik 1 1 < +00. 

Examples of all these types of motion, except the oscillatory ones, were already known by Chazy. 
The first to prove the existence of oscillatory motions was Sitnikov in |Sit60j . In his paper he 
considered the restricted spatial three body problem with mass ratio ;U = 1/2 and the three 
bodies disposed in a certain symmetric configuration, called now the Sitnikov example. Later, 
Moser |Mos73j gave a new proof. His approach to prove the existence of such motions was to 
consider the invariant manifolds of infinity and to prove that they intersect transversally. Then, 
he established the existence of symbolic dynamics close to these invariant manifolds which lead 
to the existence of oscillatory motions. 

In the planar setting, the first result was by Simo and Llibre |LS80a) . Following the same ap- 
proach as in |Mos73j . they proved the existence of oscillatory motions for the RPC3BP for small 
enough values of ^. One of the main ingredients of their proof, as in |Mos73j . was the study of 
the transversality of the intersection of the invariant manifolds of infinity. For = 0, the system 
is integrable and the stable and unstable invariant manifolds of infinity coincide. Then, they ap- 
plied the classical Poincare-Melnikov Theory |Poi901 [Mel63j to ensure that, for fj, small enough, 
the invariant manifolds do not coincide anymore and intersect transversally. Nevertheless, to be 
able to compute the Melnikov function, they considered the Jacobi constant (see ([2])) large 
enough. For these values of the Jacobi constant, the Melnikov function is exponentially small 
with respect to and therefore they were only able to prove the transversality of the invariant 
manifolds provided ^ was exponentially small with respect to the Jacobi constant. This allowed 
them to prove the existence of oscillatory motions for small enough /i. Note that the orbits that 
they obtained have large Jacobi constant, which implies that they are far from the primaries 
and, therefore, they are far from collisions. Their result was extended by Xia |Xia92j using the 



lim sup llg' 



+00 and liminf||(7|| < +00. 

t—>±oo 



t— >±CXD 
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real-analyticity of the invariant manifolds of infinity. He showed that this invariant manifolds 
intersect transversally for any mass ratio € (0, 1/2] except for a finite number of values, and 
thus he proved the existence of oscillatory motions for any mass ratio except for these values 
(see also |Mos73j for a similar argument for the Sitnikov problem). 

Following the same approach, in |MP94j some formal computations are done to show that, 
for the restricted planar elliptic three body problem, the invariant manifolds of infinity intersect 
transversally for arbitrarily small mass ratio. Taking the eccentricity small, a rigorous compu- 
tation of the Melnikov function is done in |DKdlRS12j . The existence of oscillatory motions 
has also been proven for the (non necessarily restricted) collinear three body problem |LS80bj . 
All the mentioned works follow the approach initially developed in |Mos73j . that is, they relate 
the oscillatory motions to transversal homoclinic points to infinity and symbolic dynamics. A 
completely different approach using Aubry-Mather theory and semi-infinite regions of instability 
has been recently developed in [GKIH IGKlObl IGKlOaj . This new approach has allowed the 
authors to prove the existence of orbits which initially are in the range of our Solar System and 
become oscillatory as time tends to infinity for the RPC3BP with a realistic mass ratio for the 
Sun-Jupiter pair. 

The mentioned works deal with the problem of the existence of oscillatory motions in different 
models of Celestial Mechanics. Once the existence is known, the natural question is to measure 
how abundant they are. The only result in this direction is the recent paper |GK12j . in which the 
authors study the Hausdorff dimension of the set of oscillatory motions for the Sitnikov example 
and the RPC3BP. Using [Mos73j and |LS80aj . they prove that for both problems and a Baire 
generic subset of an open set of parameters (the eccentricity of the primaries in the Sitnikov 
example and the mass ratio and the Jacobi constant in the RPC3BP) the Hausdorff dimension 
is maximal. As pointed out to us by V. Kaloshin and A. Gorodetski, the present paper and the 
techniques developed by them in |GK12j lead to the prove of the existence of a set of maximal 
Hausdorff dimension of oscillatory motions for any value of the mass ratio and a Baire generic 
subset of an open set of Jacobi constants. 

The purpose of this paper is to improve the results of |LS80aj and |Xia92j . We prove, using 
as a perturbative parameter, that the transversality of the invariant manifolds of infinity 
holds for any value of ^ £ (0,1/2] and the Jacobi constant J big enough. Moreover, we are 
find a curve in the parameter plane (/x, J) where the stable and unstable invariant manifolds of 
infinity undergo a cubic homoclinic tangency. Note that in this setting, classical perturbative 
techniques such as Poincare-Melnikov Theory do not apply because the difference between the 
stable and unstable invariant manifolds is exponentially small with respect to J^. That is, we have 
to face what is usually called the exponentially small splitting of separatrices phenomenon. This 
phenomenon was discovered by Poincare |Poi90j while studying the non-integrability of the n- 
body problem and he called it the Fundamental problem of mechanics. It has drawn considerable 
attention in the past decades but, due to its difficulty, it has essentially only been considered in 
toy models [HMSBSl IDS921 ll'regTl IGelOOl IGUSlOj or in general systems under hypothesis that 
tipically do not apply when one wants to study problems from Celestial Mechanics |DS971 IGel971 
IBF041 IBFGSll] IGual2j . For instance, most of the results deal with Hamiltonians which are 
essentially either polynomials or trigonometric polynomials with respect to the state variables. 
The present paper is the first one which proves exponentially small splitting of separatrices in 
Celestial Mechanics and also the first one which deals with an irrational Hamiltonian without 
assuming artificial smallness conditions on the parameters of the problem. 

The main result of the present paper is the following. 
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Theorem 1.1 (Main Theorem: version 1). Fix any fj, € (0,1/2]. Then, there exists an orbit 
{q{t),p{t)) of ([T]) which is oscillatory. Namely, it satisfies 



hm sup II g 



+00 and hminf||g|| < +oo. 

t— >±oo 



t— >±oo 



As we have explained the RPC3BP has as a first integral the Jacobi constant ([2]) and we 
take it big enough. Thus, in fact we have a more precise knowledge of the oscillatory motions 
we obtain. 

Theorem 1.2 (Main Theorem: version 2). Fix any n G (0,1/2]. Then, there exists Jq > big 
enough, such that for any J > Jq there exists an orbit {qj{t),pj(t)) of ([T]) in the hypersurface 
J'{q,p,t; n) = J which is oscillatory. Namely, it satisfies 



Moreover, if we take fi* G (0, 1/2), there exists J* > big enough such that, for any J > J* and 
^ E (0, /i*], the previous claims hold. 

The main difficulty to prove these theorems is to show that the invariant manifolds of infinity 
intersect transversally. To state the corresponding result, we first need to set up some notation 
and formally define what the invariant manifolds of infinity are. First, writing the system in 
rotating coordinates it turns out to be a two degrees of freedom Hamiltonian system. At each 
energy level, infinity corresponds to a periodic orbit with two dimenional stable and unstable 
manifolds. Their intersection is studied in Section [2j First, Theorem 12.11 gives an asymptotic 
formula for their difference. Then, in Theorem 12.21 we see that, for any G (0, 1/2] and big 
enough, they intersect transversally. Theorem 12.21 also gives a formula for the area of the lobes 
between two consecutive homoclinic points in an associated Poincare map. Then, Theorems ll.il 
and 11.21 follow from Theorem 12.21 using the reasoning in |LS80a| . Furthermore, we are able to 
detect a bifurcation curve in the parameter space (/i, J) where the transversality of the invariant 
manifolds is lost since a cubic homoclinic tangency appears. This fact is stated in Theorem 12.31 

The rest of the paper, that is Sections [SHSl is devoted to prove Theorems 1 2 . H [2^ and [T3l In 
Section [31 following |LMS031 ISauOlj , we look for parameterizations of the invariant manifolds as 
graphs of generating functions in a suitable domain using the Hamilton-Jacobi equation and we 
state Theorem 13. 2| which gives the difference of these generating functions. Then, Theorems 12. II 
and 12.21 follow easily from Theorem 13.21 The prove of Theorem 12.31 is deferred to the end of 
Section EJ 

Theorem l3.2l is proved in the subsequent sections. First, in Section^ we state certain analytic 
properties of the unperturbed separatrix which are crucial to obtain solutions of the Hamilton- 
Jacobi equation. Then, in Section [5j we prove the existence of "solutions" of the Hamilton- 
Jacobi equation, which are periodic in one of their variables, in certain complex domains. The 
quotation marks refer to the fact that we do not obtain actual solutions of the Hamilton-Jacobi 
equation. Indeed, to proof exponentially small splitting of separatrices, usually one has to obtain 
parameterizations of the invariant manifolds in complex domains which reach a neighborhood of 
the singularities of the unperturbed separatrix (see for instance |BFGSll| ). This is not possible 
in the present problem since the parameterizations blow up before reaching these neighborhoods 
of the singularities. Instead, we consider the Fourier series of the generating functions (see 
Section \5\ for more details) and we show that, even if the parameterizations blow up, their 
Fourier coefficients are well defined in the corresponding complex domains. This turns out to 



limsup \\qj 



+00 and liminf HgjH < +oo. 



t— >-±oo 
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be sufficient to study the exponentially small splitting of separatrices. Finally, in Section [6l we 
complete the proof of Theorem 13. 2i Analyzing the difference between two solutions (as formal 
Fourier series) of the Hamilton-Jacobi equation, which correspond to the stable and unstable 
manifolds, in the complex domains, we deduce exponentially small estimates for the difference 
between the generating functions for real values of the variables. At the end of this section we 
also prove Theorem 12.31 



2 The invariant manifolds of infinity 

To study the invariant manifolds of infinity, it is more convenient to express the Hamiltonian ([T]) 
in polar coordinates. It is given by 

H{r, a, y, G, t; /i) =— + 77^ - U{r, a - t; fj,) 

1 



- o ' T 2 U{r,a- t;fi), 

where (r, a) are the polar coordinates of the configuration space and (y, G) are the symplectic 
conjugate variables. That is, y is the radial velocity (or momentum) and G is the angular 
momentum. The function U is the Newtonian potential, which is given by 

U (r, (/>; ^i) = — ■ — — + 



(r^ — 2/ir cos (j) + M^)^^^ ('^^ + 2(1 — fi)r cos (/> + (! — /u)^)^^^ 

and therefore U satisfies U = 0{fi). Nevertheless, recall that we are considering any /x G (0, 1/2] 
and therefore n is non necessarily small. 
The associated equations are 

r = y 

G = daU (r, a — t; fi). 

Call ^t,to = (*^r,io' *^Mo' *^Mo' ^Mo) associated to this equation. Then, the stable and 

unstable manifolds of infinity are defined as 

W'^ = \{r,y,a,G) eR^ xTxR: lim '^^^^ (r, y,a,G; ^l) = 00, lim (r, y, a, G;fi)=0 
= I (r, y, a, G) G m2 ^ T X M : ^ lim (r, y,a,G; ^l) = 00 lim «>^,,, (r, y,a,G; fi) = 0\ , 



t—>-—co ' ^— 00 



(5) 



where T = M/(27rZ). It is known that these invariant manifolds are analytic (see jMcG73j ). 

It is well known that the RPC3BP possesses a symmetry which, in polar coordinates, means 
that the system Q is reversible with respect to the involution 

n{r,y,a,G) = {r,-y,-a,G). (6) 
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We will use this fact to obtain symmetric properties of the parameterizations of the invariant 
manifolds. 

When ^ = 0, the RPC3BP is reduced to a central force equation. This system is autonomous 
and therefore has conservation of energy H. Moreover, the angular momentum G is also con- 
served and then the system is integrable. This implies that the invariant manifolds of infinity 
coincide along a homoclinic manifold, = W^. This manifold is formed by a family of 

homoclinic orbits to infinity which perform Keplerian parabolic orbits. Recall that by "infin- 
ity" we mean (r, y) = (+00, 0) and then, "infinity" is foliated by periodic orbits which can be 
parameterized by the angular momentum G and are of the form 

^Go = {{r,a,y,G) : r = 00, y = 0, a G T, G = Go} . 

All the associated homoclinic orbits belong to the level of energy H = and can be parameterized 
by the angular momentum and its initial condition in the angular variable. We denote these 
orbits as 

Zh{u;Go,ao) = {ri,{u;Go),ao + ah{u;Go),yh{u;Go),Go), (7) 

and we fix the origin of time such that yh(0;Go) = and ah(0;Go) = (see Figure [1]), which 
makes the homoclinic orbit with oq = symmetric, that is, 

TZ{z^{u;Go,0)) = z^{-u;Go,0). 

Commonly, in the study of this problem, one considers McGehee coordinates |McG73j r = 
2x~'^ which send infinity to zero. In these new coordinates the system is still Hamiltonian 
with a non canonical symplectic form and the origin becomes a parabolic periodic orbit with 
stable and unstable invariant manifolds. Since the symplectic form is non canonical in McGehee 
coordinates, their use is rather cumbersome. In the present work we prefer to stick to the original 
variables, in which the invariant manifolds can be characterized as above. 

y 




Figure 1: Projection on the (r, y) plane of the separatrix d?]) of system (jH) with /i = 0. It also 
correspond to the separatrix of the Poincare map Vgo^o (flSj) . 

Before studying these invariant manifolds, we perform some rescalings that make the pertur- 
bative character more apparent for big enough angular momentum Gq. In these new variables, 
the homoclinic ([7]) will be independent of Gq, which will facilitate the exposition of the results. 
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2.1 The RPC3BP as a nearly integrable Hamiltonian System with two time 
scales 



It is well known that the existence of exponentially small phenomena usually arises when the 
system possesses two different time scales. More concretely, the exponentially small splitting of 
separatrices often appears when the system has combined elliptic and hyperbolic (or parabolic) 
behavior in such a way that the elliptic one is much faster than the hyperbolic (or parabolic) 
one. Even if this is not obvious looking at Hamiltonian this is also the case in this problem. 
To make this fact more apparent we perform a change of variables which is simply a rescaling. 
It turns out that the main theorems of this section, Theorems 12.11 12.21 and 12.31 are significantly 
simpler in the rescaled variables. 

Before showing the rescalings, let us give some heuristic ideas why this system has two time 
scales. Taking arbitrarily high Jacobi constant implies that the masless body is extremely far 
from the primaries, that is, in a neighborhood of infinity in the configuration space. Moreover, 
recall that the oscillatory motions lie close to the invariant manifolds of infinity which are 
formed by orbits whose speed tends to zero as its position tends to infinity. These two facts 
together imply that the masless body has an extremely slow motion compared with the primaries. 
Therefore, after the suitable rescaling one has a masless body which has a speed of order one, 
which is perturbed by the motion of the primaries that now are rapidly rotating. 

Let us recall that, even for system (jj]) neither H nor G are preserved as happened for the 
case /i = 0, there is still a conserved quantity: the Jacobi constant. In polar coordinates, it is 
given by 

J{r, a, y, G, t; fi) = H{r, a, y, G, t; fi) - G. (8) 

We first fix it to be ^7 = Gq. Then, in this invariant hypersurface we perform the following 
changes of variables 

r = GqT, y = G^^y, a = a and G = GqG (9) 
t = Gls, (10) 



and we rescale time as 
obtaining the system 



d_ 

Ts 

G' 1 



—r = y 
as 



d_~ _ G_ 

ds 'P 

^^G = GldaU{Gl?,5-Gls;fi). 

Now the two time scales become clear. Indeed, in these variables we have that dy/ds ~ dr/ds ~ 
1, which are the variables that will define the separatrix, whereas the perturbation dependence 
on time is fast. Calling 

V{¥,(t>;fi,Go) = GlU{Gl¥,4>;f^), 

we have that 

nr, <P; Go) = '-^ + ^ - I, (11) 



■P - 2(^)f COS + (^)2 j (P + 2(^)f cos (/> + (^)2 
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and we obtain the system we are going to study 



—r = y 
as 

d ^ & 1 _ _ 

—y = — - ^ + d^Vir, a - Gqs; /i, Go) 

as r r /-^2) 

±~ _G_ 

ds 

^G = daV{?,a-Gls;fi,Go). 

Note that V ~ fJ-G^'^ and thus, since we are taking Go ^ 1, we are dealing with a fast oscillating 
small perturbation. When one studies the sphtting of separatrices phenomena in the resonances 
of nearly integrable Hamiltonian systems typically the perturbation has the same size as the 
integrable unperturbed system (see |Tre971 IBFGSll] IGual2j ) . This is not the case in the present 
problem. This fact will facilitate the study of the difference between the invariant manifolds 
since we will not need to use inner equations and complex matching techniques or continuous 
averaging as were used in those papers. 

The rescaling ([SJ is conformally symplectic and therefore the new system p2|) is Hamiltonian 
with respect to 

„ y^G^l 

H{r, a, y,G,s;n,Go) = Y + ^-~ - V{r, a - Gqs; /x, Gq). (13) 



r 



The Jacobi constant is now J = Gq — GqG and the periodic orbit at infinity is given by 
(r, a, y, G) = (oo, 5, 0, 1), which belongs to the surface of Jacobi constant 

J {gI?, 5, Go ^y, GoG, Gls- /i) = -Gq. 

One of the main advantages of this new set of coordinates is that the parameterization of the 
separatrix of the unperturbed system ([7]) in the rescaled variables, 

(f,a,y,G) = (ni(s),ah(s),yh(s),Gh(s)), (14) 

is independent of Gq (and also of fi). In particular, Gh(s) = 1- Moreover, note that after the 
rescaling we still have 

yh(0) = and 5h(0) = 0, (15) 

and therefore 

ni(s) = rh(-s), 2/h(s) = -y\i{-s) and a^{s) = -5h(-s). (16) 



2.2 Main result: intersection of the invariant manifolds 

Let the invariant manifolds and be the rescaling of the invariant manifolds and 
in ([5]) . Their coincidence along the homoclinic manifold (jl4p for ^ = is due to the integrability 
of system ()12p for this value of /i. When /_f > 0, these manifolds do not longer coincide but they 
intersect. To study the transversality of this intersection is the main goal of this section. Before 
stating Theorems 12 . 1 1, YI72\ and [2^ we first need to explain how we measure the transversality of 
the invariant manifolds. 
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Note that the dependence on s in ()13p is only through a — Gq^s. Thus, one can ehminate 
the time dependence in the Hamiltonian by defining a new angle (f) = a — GQ^s. Then, we obtain 
the new Hamiltonian 

n{?, cP, y, G; /i, Go) = ^ - G|]G + ^ - i - V{?, <p- /x. Go), (17) 

which is an autonomous two degrees of freedom Hamiltonian and therefore has conserved energy 
H. The energy conservation corresponds to the conservation of the Jacobi constant in the original 
coordinates. In fact, ^(r, a — GqS, y, G; n, Go) = G'^J'{GQr, a, G^^y, GqG, Gqs; h). This new set 
of coordinates is usually called (rescaled) rotating polar coordinates since they are set in a 
rotating frame with respect to the primaries. Namely, in these new coordinates the primaries 
remain fixed at the horizontal axis while the massless body revolves around them. 

In these new coordinates, the invariant manifolds and are three dimensional. If 
one fixes a level of energy 'H(r, 0, y, G; Gq) = — Gq, the corresponding invariant manifolds 

W^,Go = n {nr, <p, y, G; ^, Go) = -G^} 
W^,Go = n {^(f, <P, y, G; ^, Go) = -Gj]} 

are two dimensional. Moreover, if one takes the energy high enough, in a neighborhood of these 
invariant manifolds one has that 4> ^ 0. This implies that the flow associated to the Hamiltonian 
([TT]) restricted to a level of energy '^(r, cp, y, G; n, Go) = — Gg induces a Poincare map 

'Pgo,4>o :{<A = M — ^ {0 = 00 + 27r} ^^^^ 
(r,y) ^ VGo,4,oir,y)- 

This Poincare map is two dimensional and is area preserving since the flow is Hamiltonian (more 
precisely the Poincare map preserves the symplectic form Q = dr A dy). Now the invariant 
manifolds become invariant curves 7**'^ (see Figure [2]). We focus on the part of these invariant 
manifolds with y > and we consider a parametrization of 7"''' of the form 

r = rii{v) 

y = {v;fJ',Go), 

where rtiv) is the r component of the separatrix parameterization ()14p . Then, following [SauOlj . 
to measure the distance between the invariant manifolds along a section r = const, it suffices 
to measure the difference between the functions ■ Note that the used parameterization is 
equivalent to write the curves as graphs y = y"'^(r). Nevertheless, we use the auxiliary parameter 
V since it simplifies the formulas in Theorem [2Tj We also give, in Theorem [221 a measure of the 
area of the lobes that are formed between the invariant curves since it is a symplectic invariant 
(see Figure Proceeding analogously, one could easily give an asymptotic formula for the 
angle between the curves at homoclinic points. Indeed, in Theorem 12.31 we show that if one 
take values of the parameters in certain curve r/ in the plane (;U,Go), the invariant curves 7"'* 
in , besides a transversal intersection, have a (cubic) homoclinic tangency. 
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y 




Figure 2: Stable and unstable invariant manifolds of infinity for the Poincare map "Pgc^o ™ (|18p . 

Theorem 2.1. Consider the invariant manifolds of infinity and of the system 

associated to Hamiltonian p!7j) in the level of energy T-l{r,(p,y,G; n,Go) = —Gq and the corre- 
sponding invariant curves 7"'* of the Poincare map "Pgoj^o- Then, there exists Gq > such that 
for any Gq > Gq and fi G (0, 1/2], 

• the curves 7"''^ have a parameterization of the form ()19p and 

• if we fix a section r = r* , the distance d between these curves along this section is given by 

l-2//^3/2„-gi 



d=yh{v*) nil - fJ.)y/7r 



—^G^e-^ sin (<^o - 5h(^;*) + Ggz;*) 



+ 8G^/^e-^ sin 2 (0o - a^W) + G|]r;*) 
+ O f (1 - 2^)Goe-^ + Gi]e- 



where v* is the only v > such that r{v*) = r* and y\i{v) and ah(f) are the y and a 
components of the unperturbed homoclinic ()14p . 

This theorem is proven in Section [3l Observe that the distance is exponentially small with 
respect to Go, which is taken as a large parameter. Moreover, note that the size of the first order 
in the formula of Theorem 12.11 is significantly different depending whether ^ 7^ 1/2 or /i = 1/2. 
The reason is that Hamiltonian (|17p is periodic with respect to (j) with period 2ir, in the first 
case, and vr, in the second one. The physical explanation of the different periodicity goes as 
follows. For fi ^ 1/2 (see left picture of Figure [3|) the primaries rotate around the center of mass 
with the same period 27r performing circles of different radius. On the contrary, when fi = 1/2 
both bodies have the same mass and therefore, they move along the same circle where the two 
bodies are placed in diametrally oposed points (see right picture of Figure [3|). Therefore, in 
the case = 1/2 the period of the system is the half of the period of the primaries around the 
center of mass. As usually happens in exponentially small splitting of separatrices phenomena, 
the smaller the period of the fast perturbation, the smaller the distance between the manifolds 
(see for instance |Nei84j ) . 
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nil — 1 — fi 



Figure 3: The motion of the primary bodies in the cases /.f 7^ 1/2 and fi = 1/2 respectively. 

Looking at the formula given in Theorem 12. H the zeros of the distance are given, up to first 
order, by the zeros of the function 

f{x) = (1 — 2fi) sinx + 16\/2Goe~~3~ sin23; where x = (po — ah{v) + GqV. 

The number of zeros for x G [0, 27r) and their nondegeneracy depends strongly on the relation 
between the parameters fi and Gq. If we fix / 1/2, and we take Go > big enough, f{x) ~ (1 — 
2fi) sinx and therefore we will have two nondegenerate zeros, which give rise to two transversal 
homoclinic points. On the contrary, for = 1/2 and Gq big enoug h, f{x) - 16\/2Gge-Go/3sin2a; 
and we will have four nondegenerate zeros. These facts are stated in Theorem 12.21 Clearly, 
between these two regimes the system undergoes a global bifurcation where one of the two 
transversal intersections becomes a cubic tangency where two new homoclinic points are born. 
This occurs in a curve r/ in the parameter plane (see Figure H]) as is stated in Theorem 12.31 



Go 













Go 





Figure 4: Bifurcation curve r] in the parameter space where the homoclinic tangency is under- 
gone. 

Theorem 2.2. Fix fi G (0, 1/2]. Then there exists G* > such that for any Gq > G* , 

• the invariant curves 7"'* of the Poincare map ^^0,00 'i-ntersect transver sally and 

• the area of the lobes between the invariant curves 7"'* between two transversal consecutive 
homoclinic points is given by 

(i + o(Go-/^)). 



A = fl{l - fi)VTT 



1-2/^^-3/2 



3 



+ 8G, 



2Gl 
3 
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Theorem 2.3. Let Gq > be the constant introduced in Theorem \2.1[ Then, there exists a 
curve rj in the parameter region 

of the form 

= ^r(Go) = \- IQ^Gle-^ (l + O (Go , 
such that, for {fi, Gq) E rj, 

• the invariant curves 7"'^ of the Poincare map "Pgc^o have a cubic homoclinic tangency 
and a transversal homoclinic point and 

• the area of the lobes between the invariant curves 7"'* between the homoclinic tangency 
and a consecutive transversal homoclinic point is given by 

Theorem 12.21 is proven in Section [31 The proof of Theorem 12.31 is deferred to the end of 
Section [6i 



X (GS,+oo), 



2.3 Prom transversal homoclinic points to oscillatory motions 

The existence of oscihatory motions in the RPC3BP was obtained by Llibre and Simo in [LSSOaj 
for large values of Go and exponentially small values of the mass ratio fi, concretely for fj, = 
0(^g-Go /3^^ Their arguments follow the ones developed by Moser in the so called Sitnikov 
problem in |Mos73j . The proof requires, as a first step, to control the local behavior near infinity. 
This is done using McGehee coordinates [McG73j . in which infinity becomes a parabolic critical 
point at the origin, and then the so-called Shilnikov coordinates to study the local map near 
zero 

The local behavior of the system at infinity turns out to be the same for all values of the 
parameter n G [0,1/2] and Gq big enough. The reason is that when one considers system (jH), 
makes the change of variables (j) = a — t to express it in rotating coordinates, performs the 
McGehee change of coordinates r = 2x~'^, uses the energy reduction and reparameterizes time 
s = s{0) to have dcj)/d6 = 1, one obtains (see [LSSOaj ) the system 

dx x^y x'^y 

f^=--K^- + m^-X{e;,)x^ + 0,,ix,y) 

where K = J — /i(l — ^) and \{6; ji) = ;^^(1 — At)(l — 3cos^ 9). Since the dependence on ^ 
only appears in the higher order terms, the local behavior is the same regardless the value of 
the parameter S [0, 1/2]. Therefore, the techniques developed by Moser in jMos73j that were 
applied in [LSSOaj for the case fi small enough, are also valid in our setting. 

The second step of the proof is to show that the invariant stable and unstable manifolds of 
infinity intersect transversally. It was in this step that the smallness of the parameter /i was 
needed in [LSSOaj . Theorem 12.21 guarantees this transversality for all values of /x. 
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Hence combining this fact with the local results of |LS80aj we obtain Theorem 6.1 of |LS80aj 
for any value of /U € (0,1/2], which shows that there exists a Cantor set in the phase space 
where certain return map associated to system is conjugated with a shift of infinitely many 
symbols. From this result, Theorems 11.11 and ll. 21 follow directly. 

Another way of of proving Theorem 11.11 and ll. 21 would be to consider the Poincare map ()18p in 
McGehee coordinates, which has a parabolic fixed point at the origin with transversal homoclinic 
points. Then, as Moser mentions in |Mos73] for the Sitnikov problem, one could adapt the proof 
of the classical Birkhoff-Smale theorem to parabolic fixed points with transversal homoclinic 
intersections to obtain a Cantor set of the phase space with dynamics topologically conjugated 
to a shift of two symbols. Using this idea and proving a C^-version of the Lambda Lemma, 
Gorodetski and Kaloshin in the recent preprint |GK12j . prove that the Hausdorff dimension of 
the set of initial conditions which lead to oscillatory motions is maximal. Their proof also needs 
the same requirements about the parameters /x and Go of |LS80aj because their reasoning uses 
the transversality of the invariant manifolds proved in |LS80aj . Private conversations with the 
authors indicate that, using Theorem \2.2\ their results can be extended for a larger set of the 
parameters. 

3 Parameterizations of the invariant manifolds as generating 
functions 

As we have explained in Section 12. H to prove Theorems 12. H 12.21 and 12.31 we consider the rescaled 
system in rotating coordinates given by Hamiltonian (|17|) . The main two advantages of this 
choice are the following. First, the unperturbed separatrix is independent of the parameter Gq 
and, second, working in rotating coordinates our system has two degrees of freedom (and thus 
conserved energy) instead of two and a half as happens in the original variables. 

To study the difference between the invariant manifolds of infinity, we follow the approach 
proposed in jLMSOSllSauOlj . That is, we take advantage of the fact that the invariant manifolds 
are Lagrangian and therefore they can be locally parameterized as graphs of generating functions 
which are solutions of the so-called Hamilton-Jacobi equation. 

We consider the Hamilton-Jacobi equation associated to Hamiltonian ()17p and we look for 
functions S(r, (p; fi, Gq) such that 

{y,G) = {9?^S{r,(l);fi,Go),d^S{r,(l);fi,Go)) 

parameterize the invariant manifolds as a graph. Then, the Hamilton-Jacobi equation reads 

n{?, <A, drS, d^S; fi, Go) = -Gl (20) 

Recall that we put — Gq in the right hand side since it is the level of energy where we are looking 
for the invariant manifolds. 

For the unperturbed Hamiltonian, that is, considering V = in (|17p . this equation simply 
reads 

i {(kSf - Gld^S + ^ (d^Sf -1 = -Gl 
It has a solution of the form 

5o(f,<^) = + /(?), (21) 

where / is any solution of 
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In |LMS031ISau01j . the authors deal with unperturbed separatrices which can be written globally 
as graphs. This is not possible in the present problem, as can be clearly seen in Figure [TJ since 
the separatrix has a turning point at (r, y) = (1/2,0). Thus, we will deal with this equation 
for r 7^ 1/2. This fact will carry some technical problems while proving the existence of the 
invariant manifolds in certain domains later in Section [5j 

We look for solutions of (j20p close to (j2ip . To this end, we write S = Sq + Si and then the 
equation for becomes 

d^fdrSi + ^ {d^Sif - Gld^Si + i^d^Si + ^ {d^Sif - V{?, ^u, Go) = 0. 

To look for solutions of this equation we proceed as in |LMS031 ISauOlj and we reparameterize 
the variables (r, (/)) through the unperturbed separatrix (|14p . Namely, we consider the change 

i?A) = irhiv),C + Mv))- (22) 
We define the new generating function 

r(u,^;^,Go) = 5(fh(w),C + 5h(t^);/-i,Go), (23) 
which can be correspondingly written as T = Tq + Ti where 

and 

Ti{v,^;fi,Go) = Si{rh{v),^ + 5h(u); /U, Go), 
whose associated Hamilton-Jacobi equation is 

2 



a.Ti - Gld^Ti + ^ (d,Ti - i^d^Ti) + (s^ro^ - v (fh, e + 5h; ^u, Go) = 0. 

2y^ V / 2r^ 



(24) 



Note that the change of variables ()22p implies that we are looking for parameterizations of the 
stable and unstable invariant manifolds of the form 

r = rh{v) 

y = Vhiv) + yh{v)~^ {dyT^'^v, fi, Go) - ri,{v)~^d^T^'\v, Z^, Go)) 

<P = ? + ah(v) 

G=l + ^^T^^\v,i■,^Ji,G^), 

where T"'* are solutions of equation (|24p with asymptotic boundary conditions for the unstable 
manifold 

hm y-\v)d,,T^{v,in^,G^) = Q 
lim ^^T^{v,i■,^i,GQ) = Q, 

and analogous ones for the stable manifold taking v — t- +oo. By its definition in (|lip . V has the 
symmetry property 

V{r,-e) = V{r,e). (26) 
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Taking into account this fact and ()16p . one can easily see that if Ti{v,S,) is a solution of ()24p . 
— Ti(— f, — ^) is also a solution. Thus, since T" has to satisfy the just mentioned asymptotic 
boundary conditions and the opposite ones, the generating functions parameterizing the 
invariant manifolds must satisfy 

Tf{v,0 = -T^(.-v,-0. (27) 

This means that if one is able to prove the existence of, for instance, the unstable invariant 
manifold, the existence of the stable one is guaranteed by the symmetry. 

Theorems 12.11 and 12.21 will be deduced from the study of the difference of the generating 
functions T" and Tf and from the difference between their derivatives given in Theorem 13.21 
below. The proof of Theorem 12.31 is postponed to the end of Section [6l since we will need a 
slight modification of these generating functions to deduce the existence of tangencies. 

It is worth remarking that, in view of ()15p and the second equation of ()25p . the parameteri- 
zations are undefined at = and therefore it seems impossible to have the parameterizations of 
both the stable and unstable manifold defined in a common domain. Nevertheless, in Section [5l 
we overcome this problem by using auxiliary parameterizations of the unstable manifold which 
allow us later on to recover parameterizations of both invariant manifolds of the form (|25p in a 
common compact domain contained in {w > 0} x T, which in the original variables corresponds 
to y > (see Figure [2]) . 

Before stating Theorem 13.21 we define the function that will give the first asymptotic order 
between the generating functions T"'*. This function is closely related to the so-called Melnikov 
potential or Poincare function. We define 

/ + 00 
V(ri,{v + s),C-Gls + ai,{v + s); /i, Go)ds, (28) 
-oo 

where V is the rescaled perturbed Newtonian potential defined in (jlip . Note that the classical 
Melnikov potential is the first order in fi of this function, that is 

/ di,V{ri,{v + s), ^ - GqS + ah{v + s);fi, Go)]^^^ ds. 

J — oo 

This was the function considered as a first order of the difference between the invariant manifolds 
of infinity in jLSSOaj . since they were considering perturbative methods in fj,. On the contrary, 
in the present paper fj, can take any value /U € (0, 1/2] and therefore we deal with the modified 
Melnikov potential (|28p . The following proposition, whose proof is deferred to Appendix IX\ 
gives estimates for this function. 

Proposition 3.1. The function L{v,S^; fi,Go) satisfies 

+ 00 

L{v, e; /X, Go) = (/i. Go) + 2 (^, Gq) cos £ (^ + G^v) , 

i=i 

where 

Go) = - ^(1 - ^) V^^Go^^/^e-^ {l + O {G,')) 

Go) = - 2/i(l - ^)V^Gy'e-^ (l + O (Go 
LM(^,Go) = (Gr/V^), fori>3. 
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Theorem 3.2. There exist < w_ < f+, Gq > and K > {) such that, for any Gq > Gq 
and fi £ (0, 1/2], the invariant manifolds of infinity have parameterizations of the form (|25p for 
{v,0^iv-,v+)xT. 

Moreover, the corresponding generating functions satisfy 

\Tr{v,0 - Tl{v,0 - L{v,0 -E\< K^? (1 - 2^l) Gfe-^ + KG~'^^ iJ.^e^ , 
for a constant G which might depend on ji and Gq? ^^c? 

for < m + n < 2, < m,n, where we have omitted the dependence on fj, and Gq of T^'^ , L 
and E to simplify the notation. 



3.1 Proof of Theorems [211] and [211 

We use the estimates obtained in Proposition 13.11 and Theorem 13.21 to prove Theorems 12.11 
and 12.21 The first step is to deduce parameterizations (|29p of the stable and unstable invariant 
curves of infinity 7"'* of the Poincare map VGo,(f>o in (|18p from the generating functions T"''' 
obtained in Theorem 13.21 Recall that, on the one hand, the generating functions T"'* give us 
a parameterization of the invariant manifolds of infinity associated to Hamiltonian (jl7p of the 
form ()25p . On the other hand, the Poincare section we are considering is (p = (j)Q, which has 
natural coordinates (r, y) (see (|18p ) . This implies that to obtain the parameterization ()19p of 
the invariant curves 7**'^, one has to impose ^ + 5h(w) = 0o obtaining 

+ y^{v)-' {d,T^'\v,C;fi,Go)-Uvr'd^T^''{v,C;fi,Go)) ^^^^ 



C=<^0-ah(i') 



Using the estimates of Proposition 13.11 and Theorem 13. 2[ one has that 



Yliv; /X, Go) - Yl{v;t,, Gq) =yh W {d,Tl{v,^; ^u, Go) - ^(t;,^; /x, Go)) 



5=0O-ah(") 



+ j^/iGo^/^e-^j 
=yh{v)~^d^L {v, 00 - 5h(u); Gq) 

+ O - 2^)Goe-^ + G^^V^e^ j . 

Now taking into account Proposition 13. H we obtain the formula of the distance in Theorem 12.11 
To prove Theorem [221 we first fix /x G (0,1/2). Then, taking Go > G*, G* = G*(^) big 
enough, 

Yliv; /i, Go) - Yliv; Go) = - yh(^)~V(l - /u) V^^^G^ e"^ sin (,^0 - a^^iv) + G^v) 

+ o(fiGoe-^ 
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This gives the formula of the distance for fixed fi £ (0, 1/2). Proceeding analogously, one obtains 
the following formula for fi = 1/2 and Gq big enough, 

Yliv; 1/2, Go) - Y^iv; 1/2, Go) =yh(^)~' V^G^^'e"- ^ sin 2 (<Ao - ^v) + Glv) 

i) 



We prove the transversality of the invariant manifolds and compute the area of the lobes for the 
case 7^ 1/2. The homoclinic points of infinity are just given by the equation 

Yliv;fi,Go)-Yliv;fi,Go) = 0. 

To locate the zeros of this equation we use the just obtained first order of the difference of the 
function Y^^^^{v). The zeros of this first order correspond to the zeros of sin (</)o — 5h(v) + Gqu) 
and thus are f ^ G R such that 

(30) 



bo - ah(vfc) + Gg-Ufc = fcvr, k 



which are Gq ^-close to each other. These zeros are non degenerate since for Go big enough 
5'^{v) + Gl^0. 

Using the estimates of Proposition 13.11 and Theorem 13.21 and taking Go large enough, one can 
apply the Implicit Function Theorem to show that the function Y^^{v)—Y^^{v) has nondegenerate 
zeros G M which satisfy 



^, = vl + 0{G, 



(31) 



and therefore the curves 7" and 7* intersect transversally at the corresponding homoclinic points. 

To compute the asymptotic formula for the area of the lobes, we proceed as in [BFGSll] 
and we consider the change of coordinates 

j r = rh{v) 
\y = Vhivy^w. 

Since dyrh{v) = yh{v), we have that dr A dy = dv A dw. We compute the area of the lobes in 
this new system of coordinates. Now the invariant curves are parameterized as a graph as 

w = yh{v)Y^f{v;n,Go). 

Consider two consecutive zeros which we denote by vl and vl_^_^. Then, using the definition of 
the function Y^^^ in ([29]) and that a'f^iv) = l/r^{v) (see equation (fT2]) with ^ = 0), the area of 
the lobes is given by 



A 



'k+l 



Uv) {Yl{v;fi,Go)-Yl{v;fi,Go)) dv 
{d^Tfiv, /i, Go) - ^^iv^d^T^iv, ^; Go)) 

- {d,Tl{v,C;^L,Go)-r^iv)-^d^T^iv,^;f,,Go)) , 

(rf(t;^+i, (/>o - 5h(v^+i); Go) - ^{'(4+1, - 5h(ui.+i); ^i, Go)) 
- (^i(^^fc>o - 5h(vfc);^, Go) - T'^(v\,4)Q - ai,{vl); fi,Go)) 



dv 
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Then, applying Theorem 13.21 we have that 

( _ 

A = L (wfc+i, (^0 - ahl^'fc+i); /^i Gq) - L {vl, (po - ai,{vl);fi, Gq) + O i fi^G^ 3 J . 

To obtain the formula for the area of the lobes in Theorem l2.2l it is enough to take into account 
the estimates in Proposition 13.11 that the zeros satisfy ()3ip and that the zeros of the first 
order satisfy (j30]) . 

The case fi = 1/2 can be done analogously but taking into account that now the zeros of the 
first order are solutions of 

(l)o-5y,ivl) + Glvl = k^, keZ. 

4 The integrable system and the homochnic solution 



To study the existence of the solutions of the Hamilton- Jacobi equation (|24p . we need to recall 
some analytic properties of the functions fh, and a^, which are solution of 

d ^ ^ 
—r = y 
av 

d ^_ I 1 
av 

A~ - 1 

dv ' 

close to its singularities and as \v\ — )• +oo. 

Next lemma (see |LS80al IMP94| ) gives the homoclinic solution of infinity of these equations 
satisfying (|15p after a reparameterization of time. 

Lemma 4.1. Let t{v) he the unique analytic function defined by 

2 V3 

with the convention that r is real for real values of i0. Then, the rescaled homoclinic ([1 
satisfying (|15p has the following properties: 

1. r{v) = r(T(f )), where 

2. a{y) = a{T{v)), where 



a(r) = 2arctan(r) = -ilog ( ^-p^ ) . (32) 



^In fact, 
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As a consequence, we can deduce the behavior for the homochnic close to its singularities 
and as \v\ — +00. 

Corollary 4.2. The rescaled homoclinic ()14p with initial conditions (jlSp behaves as 
rh{v) ~ 3v'^^^, yh{v) ~ 2f~"^/^, and ah{v) — ao ~ ^'^^^^^ 

as |u| — )■ +00. 

Now we give the properties of the homoclinic close to its singularities. Note that r = iti if 
and only if v = ibi/3. Furthermore, 

T±ir^ {v±i/3)^/'^. (33) 
Corollary 4.3. The rescaled homoclinic (fn|) mi/i initial conditions ([15]) behaves as 



rh(i;) ~ C ( u T 2 j , yh(^^) ~ ^ f 2 j , ah(w)~±-ln(u 



where = ±2?, and 



5 Existence of the invariant manifolds 

We devote this section to obtain parameterizations of the invariant manifolds of infinity as 
solutions of equation ()24p . Introducing the linear operator 

C{h) = d^h - Gld^h (34) 

equation ([21D becomes 

^^Ti) = --^ f - ^di,T^ - [d^T^f + U{v, i + 5h). 

where 

C/(t;,0) = y(fh(^;),^;/U,Go), 

and V is the function defined in (jlip . From now on, we omit the dependence on the parameters 
11 and Gq. To look for solutions of this equation we split the potential U as U = Uq + Ui where 

Uo{v, 0) = - /i - 3(1 - /i) cos2 0)-^^^ (35) 

and 

Ui = U-Uo. (36) 
We consider two different left inverses of C given by 

g^{h){v,^)= h{v + s,^-Gls)ds (37) 
J —00 

G'{h){v,C)= / hiv + s,^-Gls)ds, 
J +00 
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whenever these expressions are well defined. Using these operators and the map 

A:{v,0^{v,C + a{v)), (38) 

we introduce the functions 

Q* = g*{UooA), * = u,s (39) 

We will see later that Qq are indeed well defined in suitable domains. Finally, we introduce the 
new unknown Q* satisfying 

T^ = Q*o + Q*. (40) 

The equation for Q* is 

£iQ*)=J^*iQ*) (41) 

where J-'* is given by 

(42) 

We devote the rest of the section to study solutions of these equations in suitable domains. First 
in Section 15.11 we give the main ideas of the proof and define the complex domains were we will 
look for the solutions of equations (j41|) . Then, in Section [5.21 we analyze the functions C/q ^ind 
U\. Finally, in Sections 15.3115.51 we obtain parameterizations of the invariant manifolds in the 
different domains specified in Section [5.11 



5.1 Description of the proof 

The classical way to study exponentially small splitting of separatrices would be in this setting 
to look for the functions and as solutions of equations (|41|) in a certain complex common 
domain D xT„, where C C is a domain which reaches a neighborhood of size ©(Gq^) of the 
singularities of the unperturbed separatrix, namely w = ±i/3 and 

= e C/(2^Z) : |Im^| < a}. 

However, to solve equations ()4ip we have to face two different problems. 

The first one is that equation ()^T]l becomes singular at = because yh(0) = (see (fT5]) ). 
To overcome this problem we proceed as was done in |BFGS11] (see also |Gual2j ). We look 
for solutions of equations (|4ip in a complex domain for the variable v, which was called in 
that paper boomerang domain due to its shape. The important features of this domain is that 
contains both an interval of the real line of width independent of Go and points at a distance 
©(Go 3) of the singularities v = ±i/3, and that it does not contain u = 0. This domain can be 
seen in Figure [5] and is defined as follows 

Dk,5 ={v ^C] |Im-y| < tan/3iRei; + 1/3 - kGq^, |Imt;| < -tan/JiRev + 1/3 - kGq^, ^^^^ 
\lmv\ > — tan /32Rev + 1/6 — 6} , 

where k G (0, 1/3), 6 G (0, 1/12) and /3i,/32 G (0,7r/2) are fixed independently of Gq. Therefore, 
this domain is non empty provided Gq > 1. 
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Figure 5: The domains D^^^s defined in (|43p . 



To obtain the parameterizations of the invariant manifolds of infinity using equations ()4ip . 
we need to impose the asymptotic conditions 

hm y^\v)d,Q^{v,0 =0, lim d^Q''{v,O = 

(44) 

lim y^Hv)dvQ'{v,0 =0, lim d^Q'{v,O=0- 

Rc-u— >+oo Rc?)— >+oo 

These conditions do not have any meaning in the domain Df^^s since this domain is bounded. 
This implies that, to prove the existence of the parameterizations of the invariant manifolds in 
the domain D^.^j one has to start with different domains were these asymptotic conditions make 
sense and, then, one has to find a way to extend them analytically to the domain D^^^. 
Thus, the first step is to look for solutions of equations (|4ip in the domains 

D'ioo = {v^ C;Re^; < -p] 
oo,p I p; 

for some p > 0, where one can impose the asymptotic conditions (j44p . This is done in Theo- 
rem 15.51 of Section 15.31 

To analytically extend the invariant manifolds to reach -0^,5 we have to face the problem 
that these parameterizations become undefined at v = 0. To overcome it, the second step is to 
first analytically extend them to the domains 

D^;,^^^^ = |u G C;|Imu| < - tan/3iRe + 1/3 - kGq^, |Imi;| > tan/32Rei; + 1/6 - (5, 

Rev>-p'] (46) 
^p',.,5 = {^eC;-t;GZ)-*'n, 



which do not contain u = (see Figure [6]). We choose p' > p so that D^i^^g n ^ ^ 0, since 
then we can perform the analytic extension procedure. Note that D^^ s C D^p^^. This is not 
true for the unstable manifold. 
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Figure 6: The domains D^^^ and -D*^ defined in (06]). 

When one tries to perform this analytical extension procedures, the second problem arises, 
which, as far as the authors know, had not appeared before in any exponentially small splitting 
of separatrices problem. It is well known that the solutions of ()4ip become very large in the 
upper and lower vertices of the domains D^'^^, since they reach 0{Gq^) neighborhoods of the 
singularities of the unperturbed separatrix. However, in the exponentially small splitting of 
separatrices problems that have been studied in the literature, using suitable weighted Banach 
spaces, one is able to prove that, even if they become big, they can be analytically extended 
to the whole domains D^'^^. This is not possible for equations ()4ip since the solutions of 
these equations blow up before reaching these points. Indeed, if one looks for the analytic 
continuation of the functions Q"'^(f,0 = S^ez their Fourier coefficients Qu]s{v) 

grow exponentially in i near the singularities v = ibi/3 and therefore one cannot obtain these 
analytic extensions. This technical difficulty can be overcome since to study the exponentially 
small splitting of separatrices phenomena one does not need to extend the functions Q^''^{v,S,) 
but it suffices to extend their Fourier coefficients Qu]s{v)- If the functions Q'^''^ are well defined 
up to neighborhoods of the singularities, these two analytic extension procedures are equivalent, 
but in the present paper they are not. This implies that the extension theorem that we consider 
in this paper, which is explained in Section [5.4) is rather different from the ones in the literature 
(see for instance [BFGSTT]). Indeed, we will consider the Fourier coefficients of Q"^'^ as sequences 
of functions and we will study their extension in a Banach space of sequences of functions. This 
Banach space will be endowed with a weighted norm that will allow us to obtain good estimates 
for each Fourier coefficient. With this procedure we are able to analytically extend the Fourier 
coefficients of the solutions of equations ()4ip to the domai ns D^ '^ ^ even if the corresponding 
Fourier series are not convergent. This is shown in Theorem 15.101 in Section [5.41 

Once in Section 15.41 is explained the analytical extension procedure of the solutions of equa- 
tion (|4ip as sequences of Fourier coefficients, it only remains to complete the analytical extension 
procedure for the unstable manifold. Indeed, as we have already pointed out D^^ s C ^ g and 
therefore, we have the Fourier coefficients of the parameterization of the stable manifold in the 
boomerang domain (|43p . where we will study the difference between the manifolds. For the 
unstable manifold we still need to reach the points in -0^,(5 that do not belong to D^^g- That 
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is, we need to analytically extend the Fourier coefficients of to the domain 



-Ck i = \ v £ C]\lm.v\ < — tan /3iRef + 1/3 — kGq |Imu| > — tan /32Ret; + 1/6 — (5 

^ (47) 
|Imu| < tan /32Reu + 1/6 + 5 >. 



(see Figure [7]) . 





<»(l/3 - 























Figure 7: The domain 15^,(5 defined in ()T7|1 . 
Indeed, it is easy to see that 

I).,5C U 5^,5. (48) 

To obtain the analytic extension of the Fourier coefficients of Q" to the domain -Dk,5 we pro- 
ceed as follows. First, we switch to the classical flow parameterization of the invariant mani- 
folds |DS92j . which is well defined for v = This parameterization is of the form 

(f,y,(/>,G) =r"(t;,0, (49) 

where F" = Fq + F^, Fo(f,^) = ij\i{v) iVhiv) , i + oi\i{v), 1) is the parameterization of the un- 
perturbed homoclinic orbit and V^iv + t,^ — Gf^t) is a solution of the flow associated to the 
Hamiltonian (jl7|) . 

We look for this parameterization in the domain 

dI°^ = |f G C;|Imw| < -tanARe?; + 1/3 - 

' ^ . (50) 

|Imt;| < tan /32Ret' + 1/6 + 5 >, 

which can be seen in Figure [8l Finally we will switch back to the original parameterization ()25p 
in the domain D^^ s = n -Dk,<5 defined in (|47p . which does not contain w = and therefore 
equation (j^T|) is well defined. Since the inclusion pSj) is satisfied, after this procedure we will 
have the Fourier coefficients of the parameterization of the unstable manifold defined in the 
whole domain -0^,(5 • 
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7-)tt 




Figure 8: The domain D^°l defined in (07]). 



This analytical extension procedure through the flow parameterization is explained in detail 
in Section [5.5.11 It has three steps 

1. First, we derive a parameterization of the form ()49p in the domain D^^ ^ by performing a 
suitable change of variables to the parameterization ()25p , which are already known in this 
domain thanks to Theorem 15.101 This is done in Theorem 15.161 

2. Once we have obtained this parameterization ()49p in the domain D'^ ^ g, we extend it, using 
the flow, to the whole domain D^°^ . This is done in Proposition 15.201 

3. Finally, using again a change of variables, the last step is to switch back to the parame- 
terization (|25p in the domain -Dk,5 = n D^^s- This is done in Proposition 15.211 

Notation. To simplify the notation, in the forthcoming sections, we will denote by K any 
constant independent of fi and Gq to state all the bounds. 



5.2 The perturbative potential 

We devote this section to obtain estimates for the functions Uq and Ui defined in (|35p and 
respectively, in the domains defined in Section [5TTJ This estimate will be used in the subsequent 
Sections 15. 3115. 51 

First we state the following technical lemma, whose proof is straightforward taking into 
account Corollarv 14.31 



Lemma 5.1. For any v £ D, where D is any of the domains ([I3]l . ()l5]l . (06]l . (jTTp . ([50]) and 

Gq big enough, the following bound is satisfied 



1 



< 



K 



G, 



1/2 
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Lemma 5.2. For any {v,9) £ D x T, where D is any of the domains ()43p . ()45p . ()46p . ()47p . 
(|5Up and Go big enough, the functions Uq and Ui satisfy 



\Uoiv,9)\<f, 



\Uiiv,9)\<f, 



K 



Gt\Uv)? 
K 



Gl\Uv)\ 



4 ■ 



Proof. The prove of the first bound is straightforward. For the second one it is enough to bound 
the remainder in Taylor's formula using Lemma l5. II and to apply Schwarz's lemma to obtain 



U{v, 6) + 1^(1 - - 3(1 - /i) cos2 0)-^^^ 



K 



□ 



Corollary 5.3. Denoting by h^^\v) the i-th Fourier coefficient of a function h{v, 0), the following 
bounds hold 



G 



We finish this section by giving estimates for 5h(f), which will be crucial to extend the 
Fourier coefficients of the parameterization of the invariant manifolds close to the singularities 
of the unperturbed separatrix. 

Lemma 5.4. Let ^ € Z. Then, for Gq > big enough the function e^^°^^^'"^ satisfies that 
1. For V £ D, where D is any of the domains (gS]), (gSD, (gZl), 

\l\/2 





< 


v + i/3 




< 


V — i/3 



\l\/2 



fori>0 



for £<0. 



2. For V G L>^'p (see (05]) ), 



Proof. The claim follows directly from the fact that, by ([32 



„i^5h(t,(r)) _ ( i - 

i + t{v) 



and (El. 



□ 
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5.3 The invariant manifolds close to infinity 

We devote this section to prove the existence of the invariant manifolds in the domains D^^p x T^, 
where D^^^p are the domains defined in (j45p . To this end we will look for solutions of equation ()4ip 
satisfying the asymptotic conditions To obtain these solutions, we set up a fixed point 

argument inverting the operator C in (j34p in suitable Banach spaces for functions defined in 
-D^tp X To-. We do it for the unstable manifold. The existence of the stable one is given by the 
symmetry property <^T 



We first define norms for functions h : D'^^ ^ — t- C as 



\\h\\i, = sup \v''h{v)\ . 
Then, for functions h : D'^ ^ x T^- — )• C, we define 



and the functional space 

^i^,p,a = {h : D'^ p X To- — ;> C : real- analytic, ||/i||,y^o- < 00} . 

It can be checked that it is a Banach space for any fixed > 0. Moreover, since equation (|4ip 
involves both Q and its derivatives d^Q and d^Q we will need to deal with Banach spaces which 
control at the same time the norms of a function and its derivatives. To this end we define the 
norm 

mu,a = \\h\\iy,a + \\dyh\\^+i^^ + Gl\\d^h\\,y+i^^ , 
and the corresponding Banach space 

^u,p,a = {h : p X C : real-analytic, [[/iJJi.,o < 00} . (51) 

Note that 

Once we have defined the suitable Banach spaces, we can state the main theorem of this section. 

Theorem 5.5. Fix constants pQ > and (Tq > 0. Then, for Gq > big enough, there exists a 
unique solution of equation ()4ip in ^s/a^pQ^oo satisfying ()44p . Moreover, it satisfies 

for a constant bo > independent of fi and Gq. 
Furthermore, if we define the function 

LUy,0= / Ui{v + s,^-Gls + ai,{v + s))ds, (52) 
J —00 

where Ui is defined in ()36p . we have that Lf G ^5/3,po,ct(); satisfies 

m^/z,., < KfiG^^ (53) 

and 

||Q'^-L^||5A.o<^^'Go-8. (54) 



27 



The first step to prove this theorem is to invert operator C in equation ()4ip to set up a fixed 
point argument for Q". To this end we consider the operator Q = ^" defined in (|37|) . Since we 
are only dealing with the proof of existence of the parameterization of the unstable manifold, 
from now on, we omit the superindex u wherever there is no risk of confusion. 

Next lemma, whose proof is analogous to the one of Lemma 5.5 of [GOSlOj . gives properties 
about how the operator Q acts on the Banach spaces Zy^p^f,. 

Lemma 5.6. The operator Q = defined in (j57|) . when considered acting on the spaces Zy^p^^ 
satisfies the following properties. 

1. For any u > 1, Q : Zy^p^^ — )• Zy-\^p^„ is well defined and linear continuous. Moreover, 
Cog = ld. 

2. If h £ 2u,p,a for some u > 1, then 

II^WL-i,p,.<^ll^ll-,- 

3. If h £ 2^u,p,a for some u >l, then d^Q{h) G Zy^p^„ and 

\\dMh)\\y,^<K\\h\\y,^. 

4- If h £ Zy^p^fj for some u >l, then d^Q{h) G 2^u,p,ct o-nd 

\\d^gm^^^<KG^^h\\y,^. 

5. From the previous three statements, one can conclude that if h £ 2^u,p,a for some u > \, 
then G{h) £ Zy^x^p^^, and 

ig{h)Wy-i,a < K\ 



We look for a fixed point of the operator 

F = goF (55) 

where Q and F are the operators defined in ([37|) and ([^2]) respectively, in the space Z^|■l,^„^^ 
defined in (|5ip . Theorem 15.51 is a straightforward consequence of the following proposition. 

Proposition 5.7. Let us fix pQ > and (Jq > 0. Then, if Gq is big enough there exists a 
constant bo > such that the operator T in ([55]) has a fixed point £ B{bQfiGQ^) C -2^5/3.po,o-o ■ 
Moreover, the function L\ defined in ([52]) satisfies ([53]) and satisfies ([M]) . 

Before proving this proposition, we state a technical lemma which gives properties of the 
functions Qq and Ui defined in ([39]) and ([35]) . 



Lemma 5.8. The functions Qo and Ui satisfy that Qq £ Zi^pf^^a-o o,^d Ui £ ^8/3,po,o-o respectively. 
Moreover, 

||t/lO^||8/3,.o <^>Go^ 

where A is the map defined in 
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Proof. For the statements referred to Qq, it is enough to use the bound for Uq in Lemma 15.21 
and Lemma |5.6[ The statements for Ui follow from Corollary 15.31 and Lemma 15.41 □ 



Proof of Proposition 5.7. First, we see that = Q o is well defined from Z^/^ p^^^^ to itself. 
Indeed if /i € Z^/'^^p^^^^^, using Lemma 15.81 and the properties of the separatrix parameterization 
given in Corollary 14.21 one can easily see that T{h) G ^8/3 py g.^. Therefore, using the last 
statement of Lemma \5M. we obtain that Q o T{h) £ ^o- 

Now we bound J-{0) = Q o -7^(0). By the definition of T in ()42p . we haye that 

H'^) = ^{dvQ^-kd^Q^ - (a^Qof + t/i o A (56) 

Thus, applying Lemma 15.81 and Corollary 14.21 we haye that J-"(0) G ^8/3 py and satisfies 
11-^(0) 1 1 8/3,0-0 < K^iGq^. Then, applying the last statement of Lemma \5M we obtain that there 
exists a constant 60 > such that 

||-F(0)||5/3,.„<|/.Go^ 

Now we haye to prove that the operator T is contractive in BibQ^G^^) C 2^5/3,^0,0-0- Take 
/ii, /i2 G B(bofiGQ^) C 2^5/3, po,o-o- Using the last statement of Lemma 15.61 we have that 

mh2) - -F(/li)|5/3,,„ < K \\T{h2) - Hhl)\\s/3,., . 

Thus, we only need to bound the right hand side of this formula. To this end, we write it as 
- 4-2 {{'^(^vQo + dvhi + (9„/i2) - ^{2d^Qo + d^hi + (9^/12) ) {d^h2 - d^hi) 



2r^ 



(57) 



Using the bounds for Qq obtained in Lemma 15.81 Corollary 14.21 and that /ii, ft-2 G B{bofiGQ ^) C 
'25/3,po,^o' ^6 obtain that 

mh2) - Hhl)\\s/3,.o ^ - ^l||5/3,.o- 

Then, since /ii, /12 G 25/3,po.o-o> conclude that 

mh2) - -F(/ii)|5/3,,„ < K^Go^||/l2 - /ll||5/3,.o < ^/"Go''i[/i2 " /ij5/3,.o. 

and therefore, taking Go big enough, it is contractive and has a unique fixed point in B{bofxGQ^) C 
-2'5/3,po,o-o- This completes the proof of the first claim. 

Since = Q{Ui o A), inequality ()53p follows from Lemma 15.81 and the second statement of 
Lemma [53 
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Finally, to prove inequality ()54p we note that, since is a fixed point of J-, 

HQ" - L?ii5/3,<xo < w^m - mh/z,., + 11-^(0) - L^'iisA.o- (58) 

The first term in the right hand side can be bounded using that J- is Lipschitz, 

To obtain a bound of the second term in the right hand side of (|58p we first note that, by 
Lemma 15.81 

11^(0) -^lO^||8/3^,^ = 

The claim follows then applying the second statement of Lemma 15.61 to -7-"(0) — L" = Q{J-{Qi) — 
UioA). □ 

5.4 The invariant manifolds close to the singularities 

As we have explained in Section 15. 1|, once one has obtained the existence of the invariant man- 
ifolds close to infinity, the usual next step when one wants to study the exponentially small 
splitting of separatrices phenomena is to extend the parameterizations of the invariant mani- 
folds to the domains D^'^g defined in ()36]l . That is, up to points at distance 0{Gq^) of the 
singularities of the separatrix (jl4p . This is not possible in this problem since the invariant 
manifolds parameterizations blow up before reaching these points. Indeed, if one looks at the 
analytic continuation of the function Lf{v,S,) = X^^g^ -^f^(^)^*'^^ defined in (j52p . which is the 
first order approximation of the function Q^{v,£,) obtained in Proposition 15.71 one can see that 

\£] 

its Fourier coefficients grow exponentially in i. Therefore, one cannot obtain this analytic 
extension and neither the one of Q". We overcome this difficulty analytically extending the 
Fourier coefficients Q'^^ instead of the function Q". It turns out that this is sufficient to study 
the exponentially small splitting of separatrices. To this end, we consider the Fourier coefficients 
of as a sequence of functions and we study their extension in a Banach space of sequence 
of functions. We endow the Banach space with a weighted norm that allows us to obtain good 
estimates for each Fourier coefficient. 

The structure of this section goes as follows. First in Section 15.4.11 we define the Banach 
spaces for the Fourier coefficients and their sequences. We also show that these spaces have 
an algebra-like structure with respect to the classical product of Fourier series. Note that the 
"Fourier series" we are dealing with are formal and therefore this algebra structure is not obvious. 
Then, in Section 15.4.21 we set up a fixed point argument and prove the existence of the Fourier 
coefficients of the parameterization of the invariant manifolds in the domains D^'^^. We deal 
only with the unstable manifold since the existence of the stable one is given by the symmetry 



1 / 1 \^ 1 

2y£ V K J 



< K^fi^Go' 



R/H.rrn 



5.4.1 Weighted Fourier norms and Banach spaces 

We devote this section to study the weighted Fourier norms. To this end we first define norms 
for the Fourier coefficients, that is, for functions h : D'^^^ — t- C. We define 

||^IU_,r.+ = sup \{v-i/3y+{v + i/3y-h{v)\ (59) 
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and the corresponding Banach space 

Xu-,y+,p,K,& = {h: Z)^,, ,5 C : analytic, ||/i|| < 00} . (60) 

Next step is to define a weighted norm for sequences of functions. To denote a sequence {/I'^^j^gz 
we keep the Fourier series notation 

but we want to stress that this series is a formal series and that h{v,S^) is not necessarily a 
function. That is, each Fourier coefficient is an analytic function defined in ^ g but its sum 
does not need to be convergent for any £T„. We define the following norm for these sequences 
of Fourier coefficient^ 

\\h\\ua = y\W'^ el^l^ (61) 

and the corresponding Banach space 

^+£/2,u-e/2,p,K,5^ \\h\\u,K,a < 00 > . (62) 



Next lemma gives algebra properties for these Banach spaces. Its proof is straightforward. 
Lemma 5.9. The spaces yu,p,K,s,iT satisfy the following properties: 

• If h £ yu,p,K,S,cT and g G yrj,p,K,S,(T, then the formal product of Fourier series kg defined as 
usual by 

{hgfHv) = Y.''^'W'^'^ 

kez 

satisfies that kg £ yu+ri,p,K,s,a- o-nd 

\\hg\\y+r,,<T < \\h\\uA9\\v,<T- 

• If he yu,p,K.,S,a, then h G yu+ri,p,K,s,a with r] > and 

\\h\Ur,,. < K\\h\\,,^. 

• If he yiy,p,K,6,a> ^^^^^ h ^ yy-r,,p,K,8,a with 7] > and 



< KG\ 



^An equivalent norm is defined by 

^Vll/iWe-"^"!! el*^!", where = sup \{v - i/S)" {v + i/S^ h{v)\ 

and 5h has been defined in (|14p . 
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In Section 15.31 we needed to control at the same time the size of functions h G ^u,p,a and 
the size of their derivatives d^h and d^h. This fact has prompted us to deal with the Banach 
spaces ()5ip . To perform the Extension Theorem 15.101 we will also need to control the size of the 
derivatives for sequences of Fourier coefficients h G yu.p.K,5,a- The derivatives of sequences are 
defined in the natural way 

Thus, we define the Banach space of formal Fourier series 
where 

= ll^l|v,<7 + \\dvh\\^+l^a + Gl\\d^h\\^+l^a- 



5.4.2 The fixed point argument 

In Theorem 15.51 we have obtained the existence of the function in the domain -D!^ x Tq-q 
and, therefore, we have obtained the existence of its Fourier coefficients Q'^^ in the domain 
^oo po ■ devote this section to obtain the analytic continuation of these Fourier coefficients 
to the domain D^_^ defined in (06]), where we choose pi > po so that the domains D^ ^^ and 
^Pi K 5 overlap. The main theorem of this section is the following. 

Theorem 5.10. Consider po and cjo the constants given by Theorem 15.51 pi > pQ, hq,6o > 
and the function obtained in Theorem \5.5[ Then, for Gq big enough, its Fourier coefficients 
can be analytically extended to D^_^ and the sequence given by the Fourier series 

satisfies that G 3^i,pi,Ko,i5o,o-o '^^'^ 

for a constant bi > independent of p and Gq . 

Moreover, the Fourier coefficients of the function Lf defined in ()52p can be analytically 
extended to D"^^ and its sequence belongs to 3^i,pi.fto,<5o,o-o- Furthermore, it satisfies 

11^1 Ikao < KpG^"" (64) 

and 

HQ" - ^5^111,-0 <^/^'G'o-^ (65) 
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We devote the rest of the section to prove this theorem. To this end, we derive a fixed point 
argument from equation (|1T|) . We follow the approach considered in [BFGSll] . We consider an 
operator G which acts on yu,p,K,s,a and is a left inverse of the operator C in (fM]) . G is defined 
acting on the Fourier coefficients. We consider vi,vi G C the vertices of the domain D^^^ (see 

Figure [6l the bar denotes the complex conjugate). Then, we define Q as 



where its Fourier coefficients are given by 

g{hp{v) = r h^°\t)dt 

J-p 



for £ < 



for e>Q. 



(66) 



Observe that the definition of the operator Q depends on the domain, since it involves the 
vertices vi, vi and also p. Next lemma, whose proof is analogous to Lemma 5.5 in [GOSlOj . 
gives properties of this operator acting on the space yu^p^K,s,a- 



Lemma 5.11. The operator Q in ()66p satisfies the following properties. 

1. If h e yu,p,K,5,a for some v > 0, then Q{h) G yu,p,K,5,cr, ^ ° G{h) = h and 

\\mh,.<K\\h\i,^. 

2. If he yu,p^K,s,u for some v > 1, then Q{h) £ yu^i,p,K,s,cT and 



Q{h) 



<K\\h\U. 



3. If he yy.p.ii^s^o- for some u >1, then dvQ{h) G yv,p,K,s,<7 and 



d,g{h) 



< K 



4- If he yu,p,K,5,a for some u>l, then d^Q{h) G yu^p^n,5,a and 

d^g{h) 



— '^ll"'l|t^,o-- 



5. From the previous three statements, one can conclude that if h £ yu,p,K,5,a for some u > 1, 
then Q{h) G yy-i,p,K,5,a and 

lQ{h)h~i,a < K\\h\\,,^. 
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We use the operator Q to set up a fixed point argument from equation (|4ip. Nevertheless, 
since we want to obtain the analytic continuation of the Fourier coefficients obtained in The- 
orem 15.51 we impose certain initial conditions. Since the operator Q involves integration from 
different initial points depending on the harmonic, we impose different initial conditions for each 
harmonic. We define 

Recall that, since pi > po, vi,vi,—pi G ^to,po ^'^d therefore, by Theorem I5.5[ F is already 
known. Moreover, it is straightforward to see that F G 3^o,pi,ko,i5o,o-o 

iFh^a, < KfiG^^ (67) 

With this prescribed initial condition, it can be seen that the solutions of equation 

Q{v,0 = F{v,0+GoF{Q){v,0, 

understood as an equation of (not necessarily convergent) Fourier series with analytic coefficients, 
give the analytic continuation of the Fourier coefficients Q^^\v) obtained in Proposition 15.71 
Note that it is not obvious that this equation is well defined for Q £ yu,pi,Ko,So,(To for any u > 0. 
Nevertheless, the algebra properties stated in Lemma \5M ensure that it is the case. 
We look for a fixed point of the operator 

S{h) =F + goF{h). 

Next proposition gives its existence and uniqueness. From it, using that the function Q of 
Theorem 15.51 is also a fixed point of S in the overlapping domain, one deduces Theorem 15.101 

Proposition 5.12. Consider po and ao the constants given by Theorem 15.51 pi > po, > 

and 6q > 0. Then, for Gq big enough, there exists a constant 6i > independent of Gq and p 
such that the operator S has a unique fixed point G B{bipGQ^) C yi,p-i^,Ko,So,'^o- 

Moreover, the (not necessarily convergent) Fourier series L" in ()52p satisfies ()64p and 
satisfies ([65]) . 



To prove this proposition we first state the following technical lemma. 

Lemma 5.13. The functions Qo and Ui, defined in (|39p and (j36p respectively, satisfy that 
,Ko,<5o,o-o '^^^ G 3^2,pi,fto,<5o,o"o- Moreover, 

lQo]\l/2,ao < KpG^^ 

||C/io^||2,<.o < A>Go^ 

where A is the map in ()38p . 

Proof. For the first bound, we can write the function Qq in (|39p as 
where 

go(«, e) = ^ e*^Go(--i)QoM(^,)e^^« + Q,y^\-p,) + ^ e^^«o(-^i)QoM(i;,)e*^«. 
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We bound the two terms. For the first one, let us point out that vi,vi, —pi G D'c^,p(, therefore 
we can bound the function Qq in these points thanks to Lemma 15.81 Then, we can deduce for 
the function Qq the bound HQoH 1/2,0-0 — KpG^^. For the second term, we use Corollary 15.31 
and Lemma [5.41 to obtain ||C^o ° ^||3/2,o-o — KpGq'^. Then, applying Lemma [5. Ill we obtain that 

For the second bound it is enough to apply Corollary 15.31 and Lemma |5.4[ □ 
Proof of Proposition \5.1^ First, we see that S is well defined from 3^i,pi,Ko,<5o,o-o *° itself. Indeed 

if he yi. 

:Pi,Ko,<5o,o-o) using Lemmas [5T9] and [5TT3] and the properties of the separatrix parameteri- 
zation given in one can easily see that T{h) G 3^2,pi,/to,'5o,o-o' where is the operator defined 
in (j42p . Therefore using the last statement of Lemma l5.11l we obtain that QoT[h) G 3^i,pi,Ko,(5o,o-o- 
Then, using also (f67|) . we can deduce that S{h) G 3^i,pi,Ko,(5o,o-o- 
We bound the first iteration 

smv,i) = F{v,i) + g o Fmv,i). 

Recall that J^(0) has been defined in (|56p . Applying Lemmas 15.91 and 15.131 and Corollary 14.31 it 
is easy to see that -F(0) G 3^2,0-0 and 1 1 -7^(0) || 2,0-0 — K^Gq^. Then, applying the last statement 
of Lemma 15.111 and taking into account (j67p , we obtain that there exists a constant 61 > such 
that 

||5(0)||i,., < l^xGo^ 

Now we prove that the operator S is contractive in BihipG^^) C 3^1,^1 , fto,<5o, 00 • Take ^i,/i2 S 
,Ko,5o,o-o ^^'^ recall that 

S{h2)-S{hi) = g{F{h2)-F{hi)). 

We start by bounding T{h2) — J-{hi). To this end, we use formula (|57p . Then, using Lemmas 15.91 
and 15.131 the properties of the separatrix parameterization given in Corollary 14.31 and that 

,pi,Ko,<5o,o"o ' w^ obtain 

\\Hh2) - -F(/il)|l2,oo ^ ^>Go'/\/l2 - /ijl,oo. 

Then, since /ii, /12 G 3^i,pi,Ko,<5o,o"o and using Lemma 15. IH we can conclude that 

and therefore, taking Gq large enough, it is contractive and has a unique fixed point in B{bifj,GQ^) C 
,pi,Ko,5o,o"o- This completes the first part of the proposition. 
To prove the second part, we introduce 

where L'^ was defined in ([52]) . By ([53|) . Fi G 3^i,pi,fto,'5o,o-o ^^"^ 

iF^t,ao<Kf,G^^. 
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Then, since their Fourier coefficients coincide, we have that 

L« = Fi + g(C/io^). 

Hence, inequahty (|64p follows from the estimate for Fi just obtained, Lemma 15.131 and the 
second statement of Lemma 15.111 

Finally, to prove inequality (|65|) we note that, since is a fixed point of S, 

HQ" - ^^Ik.o < \\S{Qn - 5(0)||i,^„ + ||5(0) - L^'lli,.^. 
The first term in the right hand side can be bounded using that S is Lipschitz, 

||5(Q«)-5(0)||i,.„<15(g")-5(0)|i,.„ 



,O"0 



To obtain a bound of the second term in the right hand side we note that 

5(0) -L'{ = F-Fi+ g(^(0) - Ui o A). 

Then, by 



and, by Lemma 15.131 

11^(0) -[/ioyl||2,.„ 



2,f70 



The claim follows then applying the second statement of Lemma 15.111 to Q{J-{Q) — Uio A). □ 

5.5 Extension of the parametrization of the unstable manifold by the flow 

From Theorem 15. 101 we have a formal parameterization of the unstable invariant manifold as not 
necessarily convergent Fourier series in the domain D^'^ . This parameterization is given by 



y 

t 



where 
with 



V 

= To + QJ^ + Q", 



(68) 



(69) 



5o was introduced in ([21]) , Qq in and is given by Theorem 15.101 Analogously for P** . 

However, since yh(0) = 0, it is no possible to extend these parameterizations to a common 
domain containing a real interval of width independent of Gq^. To compute its difference is 
necessary to have the parameteritzations of both manifolds defined in a common (real) domain. 
To overcome this difficulty we extend the unstable manifold using a different parameterization 
as has been explained 15.11 
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5.5.1 From Hamilton-Jacobi parametrizations to parametrizations invariant by the 
flow 



The first step is to look for a change of variables of the form 

Id + g : ^ {v + gi{v,0,S. + 92(.v,0), 9 = {91,92), 

in such a way that, applied to (|68|) . r** = o (Id + g) satisfies 

«>t(f"(t;,0) = f"(7; + t,e-Ggt), 

where $f is the flow associated to the Hamiltonian system ()17p . Denoting by X the vector field 
generated by this Hamiltonian, this equation is equivalent to 

£(f") = X o f (70) 

where the operator C defined in ()34p is understood to be acting on each component of F". 
In fact, writing X = Xq + Xi, with Xq the vector field for ^ = 0, and denoting by Tq the 
parameterization defined by equations (f68|) replacing by Tq, it is an immediate computation 
to see that 

£(ro) = Xo o To. (71) 

Hence, the condition one needs to impose on g to ensure that F" = o (Id + g) satisfies ()70p is 

a.F" o (Id + g){l + C{gi)) + d^f^ o (Id + g){-Gl + £(92)) = X o F" o (Id + g). (72) 

We remark that, unlike the Hamilton-Jacobi equation, this is a system of four partial differential 
equations. However, the symplectic structure implies that it is only necessary to solve two of 
them, which we choose to be the first and the third ones. Then, the other two are also fulfilled. 
By the definition of Tq in (f69|) , ([68]) and using ([7T]) , equation ([72]) is equivalent to 

C{g) = Fo{ld + g), F=(^I) (73) 

where 



and 



= Qo+ Q"- 

We emphasize that equations ([72]) or their equivalent ([73]) have to be understood at the level 
of formal Fourier series, since T" and, consequently, F**, are this kind of objects. In particular, 
the change \d + g will be a formal Fourier series. Hence, it is necessary to give a meaningful 
definition of the composition of formal Fourier series. This will be accomplished taking into 
account that the formal change of variables Id + 5 is close to the identity and using Taylor's 
formula. 

In order to find the change Id + g, we define for g = (51, 52) £ yu,p,K,5,(7 x yu,p,K,5,cr, where 
yu,p,K,&,a is the Banach space defined in ([62]) . the norm 

1 1 5 1 1 1^,0" ~ 1 1 51 Ili/jCr 92\\v,CF- 

The following two technical lemmas summarize the properties of the composition of formal 
Fourier series we need. The second is an immediate consequence of the first. 



37 



Lemma 5.14. Fix constants a' < a, 5' < 6, p' < p and k! > k such that (log k'— log k)/2 < a— a' 
and take h G yu,p,K,S,u- Us derivatives, as defined in satisfy d^^d^h G yu,p',K',5',a' o.nd 

< (^) \,,_,)„(,_°r';;l'g>^,_i„^,);,)J I'-ll..- 

Proof. Fix constants a' < a, 6' < 6, p' < p and k! > n such that (logK;' — logK)/2 < a — a' and 
take /i G 3^i,,p,K,5,<x- By dM]), 

and, by Cauchy's formula, 



{z - vY^+^ 

where 7^, is a curve with index 1 with respect to v in D"^ ^ ^. In particular, since we will measure 
the norm of this function in D"^ ^ ^ and for Gq big enough, one can choose that 5—5' > {k' — k)Gq'^ 
and p — p' > {k' — k)Gq ^ , and we can take 7^, : [0, 27r] — )• C : t 1— )• v + (k' — K)GQ^e**. Hence, 
by dSS]) 



= sup 



< 



where we have used that 



sup 



v±i/'i 



u+\i\/2 



v + {k,' -K)GQ-^e'^±i/2, 



K 

< —. 

K 



Then, by ([61]), since h G >',y,p,K,<5,CT, 



< 



ml 



u+e/2,iy-e/2 
K ^ 



< 2- 



(k' - K)'"(cr - ct' - (log k' - log k) /2)" V 1^ 
where we have used that, for 6 > 0, J2e>o'^"'^~^^ ^ 2n!/6". 



□ 



Lemma 5.15. Given g = ((71, 52), ^^t \d + g : {v, ^) 1— )• (u + gi{v, + g2iv,0)- We define the 
formal composition of formal Fourier series 

ho{U + g){v,0 = Hv + gi{v,0,C + 92{v,0) 



m=0 n=0 



Fix constants a' < a, p' < p, 5' < 5 and k' > k such that (log k' — logK)/2 < a — a' . Let 
9 = min{K' - K,a - a' - (log n' - log k) /2} > 2?? > 0. Then, 
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U h £ yu,p,K,5,a, g = (51,92) G yQ,p' ,k' ,5' ,a' X >'o,p',k',<5',<7' and WgWo^a' < vGq ^, we have that 
h = ho (Id + g) satisfies h G yv,p',K,',5',cT' o,nd 

\h\ 



Moreover, if \\g\\o^a' , ll^llo.o-' < flG^ , then f = ho (Id + g) — ho (Id + 5) satisfies 



2G3 



2r/ 



- 1 



u,a\\g — g\w,a' 



Theorem 5.16. Let pi, 5q and kq and ctq be the constants given by Theorem \5.10[ Let ai < do, 

P2 < Pi, ^1 < ^0 and Hi > HQ such that (log ki — log k,q)/2 < ao — ai be fixed. Then, for Gq big 
enough, there exists a (not necessarily convergent) Fourier series g = (91,52) G 3^o,p2,ki.5i.o-i ^ 
,p2,K,i,5i,ai satisfying 

\\g\\o,ai < ^2AiGo'^, 
where 62 > is a constant independent of /j, and Gq , such that 

f'^ = f"o(Id + 5), 

satisfies ([TOj). 

Proof. The function g is found as a solution of equation ()73p . The following auxiliary lemma 
provides the appropriate bounds on the function F introduced in ([75]). 

Lemma 5.17. F G 3^i/2,pi,«o,<5o,ao ^ ^1/2 ,P1,KO,<50,0-0 '^^^ l|-^lll/2,(T0 

Proof of Lemma 5.17 It follows directly from the bound on Qq in Lemma 15.131 and the bound 
on Q" in Theorem 15.101 □ 



Using the operator Q in ()66p . we rewrite equation (j73p as 

g = g o F o (Id + g) 



(74) 



and we look for g as a fixed point of this operator. 

Taking Gq large enough, by Lemmas 15. 171 and 15. 151 the map g ^ F o (Id + g) is well defined 
and sends the ba ll B{KiiGq^) C J^cpa.^i.^i.ai x 3^o,p2,«:i,5i,<7i to 3^o,p2ai,5i,<7i x J^o.paAi.^i.'ri • Also, 
by Lemma 15.171 



Fo (Id + g)|3=o 

which implies, by Lemma I5.11|, that 

' QoF 



< 


F 


1/2,(71 





0,0-1 



1/2,(70 



< K^lG^^. 



< KpGq\ 



It only remains to see that the map defined by (j74p is Lipschitz. But this is straightforward, 
since, by Lemma I5.15| for any g,g £ B(^K fiG^^) C 3^o,p2,«;i,5i,(7i ,P2,Kl,<5l,(7l , 

Fo(ld + g)-Fo(ld + g) < KpG^' \\g - g\\o,., , 

1/ 2,(71 

where the constant K depends on the reduction of the domain. Finally, using the last inequality 
and Lemma 15.111 we have that 



g{F o{U + g)-Fo (Id + g)) < KKpG^' \\g - g\\^ „ 

0,(71 ' 



□ 
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5.5.2 Analytic extension of the unstable manifold by the flow parameterization 

Now we can extend the flow parametrization of the unstable manifold given by Theorem 15.161 
whose Fourier coefficients are defined, up to now, for V G ^'pj.^^i,.?!' *° defined in ([50]) (see 

Figure [8]). This extension is obtained using the flow of the Hamiltonian system (|17p . 

We look for a parametrization of the unstable manifold F" satisfying (|70p . defined in -DJ?"™^ 
and which coincides with the one given by Theorem [5l6] in D'^^^^^ g^ n 0^°"^^^. 

Using the notations of the previous section, where X = Xq + Xi is the vector field generated 
by the Hamiltonian ()17p . and introducing F" such that F" = Fq + F^*, where Fq is obtained 
from ([BS]) with Tq instead of T", and satisfies C(Tq) = Xq o Fq, equation ([70]) becomes 

C(T^,)=TiT^,), (75) 

where 

£(F) = £(F) - DXo{To)T (76) 

and 

-F(F) = Xo(Fo + F) - Xo(Fo) - DXo{To)T + Xi(Fo + F). (77) 

We find a solution of this equation satisfying the appropriate initial condition. 

Like in the previous section, equation (j75p is a set of four partial differential equations. Also, 
since the initial condition is given by a formal Fourier series, we look for solutions in the space 
of formal Fourier series associated to this domain. For this reason, we introduce the space of 
formal Fourier series 

= I hiv,0 = h^'KvV^ ■■ h^'^ e {^.+e/2,u~i/2,.,5y , \\h\Ua < oo I , (78) 



where 



h : Dlj ^ C : < oo} , 



the norm || • is the one defined in (|59p . but taking the supremum in D^g" and, 

4 

\\h\\v,„ = ^ ||/ii|U,(T, 

i=l 

where the norm || • in the right hand side above is the one in (|62p . We remark that, since 
zbi/3 are at a distance 0(1) of the domain the norms || • \\u,cj are equivalent for all v. 

Therefore, we only work with v = 

Let ^{v) be any fundamental matrix of the linear system 

i{v)=DXo{Mv,i))z{v), t;Gl)^°X- 
In fact, the system above does not depend on ^. As a consequence, the matrix ^ satisfies 

= DXq{Vq)-^. 

Since Fo(f,^) is well defined and bounded for v G -0^?°™^ and -0^!°™^ is bounded, we have that 
there exists K > Q such that 



sup ui^^{\\^{v)\\q^,,\\^{vY\\<K, (79) 
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in any matrix norm. 

We choose V2 at the top of the domain -DJ|°™ and define Q hke in (j66|) . Lemma 15.111 applies 
without further modification. Then, it is a easy computation to check that a left inverse of the 
operator C in ([75]) is 

^(r) = ^g{^~^T). (80) 

Since we want to obtain an analytic continuation of the parameterization of the invariant 
manifold given by Theorem 15. 16] we introduce an initial condition defining 

+ ^(t;)^-i(-p2)rfl(-p2). 

\£] - 

Notice that the coefficients r!^ are already known at the points V2, V2 and —p2- The following 
lemma collects the properties we need about F^. 

Lemma 5.18. The function satisfies 

• £(F;) = and 

• F? G 2:o,Ki,5i,<:ri with ||F5||o,ai < KfiG^^. 

Proof. The first statement is straightforward. For the second one, it is enough to write 

= f " - Fo 

F" o (Id + 5) - To o (Id + g)) + (Fo o (Id + g) - Fq) . 



Now, using the formula of F" in ()68p . the estimates given in Theorem 15. 101 the estimates for Qq 
in Lemma 15.131 and the ones for g given in Theorem 15.161 Lemma 15.151 for the composition and 
the fact that the domain 1?^°^ only contains points at a distance of order one of the singularities 
V = =ti/3, one obtains the estimate of the second statement. □ 

We rewrite equation (j75p using the function Fj* and the operator Q in (|80p as 

F^ = F; + goJ^(F^). (81) 

It is important to remark that the definition of J- in ()77p involves the compositions Xq{Tq + T'^) 
and Xi{Tq + F") of formal Fourier series. Like in the previous section, in Lemma I5.15| this 
composition is defined through a Taylor series of the form 

X{To + F) = ^ iz)^'X(Fo)F®^ (82) 

fc>0 

However, this case is much simpler, since the vector fields Xq and Xi and the function Fq are 
true functions, analytic and appropriately bounded in the domain . 

Lemma 5.19. Let K > 0. Assume F,F G ^o,ki,<5i,o-i with ||F||o,o-i , ||F||o,o-i < KiiGq^. Then 
there exists K' > such that, if Gq is large enough, 
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• defining Y(r) = Xo{Tq + F) — Xo(ro) — DXQ{TQ)r , we have that Y(r) £ -Ho,ki,(5i,cti with 

\\Y{T)\\o,.,<K'^,G^^ 

• Xi(ro + r) G ^o,KiA,ai with ||Xi(ro + r)||o,^i < K'fiGo\ 
. ||y(r) - y(r)||o,., < A^Go^yr - r||o,.i, 

. ||Xi(ro + r)-Xi(ro + r)||o,.i < K'fiG^^\\r -f\\o,a,. 

Proof. Cauchy estimates imply that 

||Z)'=Xo(ro)||o,. < Kkl, ||Z)'=Xi(ro)||o,<x < KfiGQ^kl 
The claims follow then from formula ()82p . □ 
Now we can claim 

Proposition 5.20. Let ki, 6i and ai be the constants considered in Theorem \5.16l Then, there 
exists 63 > such that if Gq is large enough, the fixed point equation (|8ip has a unique solution 

€ B{b3fiGo^) C Zo^^.^Su.^- 

Proof. By Lemmas I5.18t 15.191 and the algebra properties of ^o,ki,5i,o-i' loaap JC : T ^ + 
G o J^(r) is well defined from B{K fiG^'^) C Zo^Ki,5i,cri to ^o,ki,(5i,(ti for any K > 0. Also, by 
Lemmas 15.18^ \57L9\ and [5.111 and using the bound (j79p on the fundamental matrix $, we have 
that 

ll^(o)llo,a, = l|r? + g{x, o ro)||o,., < ^ji^G^'. 

for certain 63 independent of Gq and /U. 

Finally, Lemmas I5.19l and l5.11l imply that IC is Lipschitz with Lipschitz constant K^Gq'^. □ 

5.5.3 From flow parameterization to Hamilton-Jacobi parameterization 

Finally we apply a change of variables to the parameterization F" = Fq + F]* obtained in Proposi- 
tion [S^Ql which is invariant by the flow, in order to obtain an extension of the parameterization 
of the form ([68]) to the domain -Dk2,(52 defined in (jl7|) for some K2 > ki and 82 > 61 (see Fig- 
ure [71). Equivalently we obtain an extension of the Fourier coefficients of T", solution of the 
Hamilton-Jacobi equation. Like in the previous two sections, since this change of variables has 
to be found in a domain which is far from the singularities v = =tz/3, the procedure is rather 
simple. 

We will find a change of variables of the form Id + /, with / = (/i, /2), such that 

7rior'o(Id + /)(t;,e)=fhW, 

7r3or'o(Id + /)(i;,0 = 5h(^;)+C, 

where vTj denotes the projection on the i-th component. We will see that this change of variables 
is unique under certain conditions. Therefore, in D^^^ 5i '^-^k2,<52 the change Id+/ is the formal 
inverse of the change ld + g obtained in Theorem 15.161 Then, 

vr2 o F" o (Id + f){v, = Mv)-' {dvT^iv, - ?i,{v)-^ d^T^ {v , 0) 
n,of^o{U + f)iv,0 = d^T^{v,0- 
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These equalities provide an analytic extension of the Fourier coefficients of T" given by Theo- 
rem EHUl 

Taking into account that vri o ro(f,^) = rh(v) and ir^ o ro(v,^) = 5h(f) + ^ and that 
= Fq + F", equations are equivalent to 

f = r{f), 

where 



-y^\v) iri,{v + Mv,0) - Mv) - Mv)fi{y,0 - vri o F^^ o (Id + f){v,0) 
Mv) - + Mv,0) - vrg o F^ o (Id + f){v,0 



Then, defining the space Z,y^f^^s,a analogously to (f75|) . with only two components / = (/i,/2) 
and the same norm, one has 

Proposition 5.21. Consider the constants ki, 5i and ai given by Proposition \ 5.20i and any 

K2 > Ki, 62 > 5i and 02 < (Xi. Then, 

• There exists 64 > independent of Gq and n such that, if Gq is large enough, the operator 
V has a unique fixed point / = (/i, f2) £ -2^0,k2,52,o-2 

ll/llo,.. < ?>4MGo^ 

• The change Id-|- / is the inverse of the restriction of the change given by Theorem \5.1b\ to 
the domain -Dp^Ai,^! ^ ^'^2,62- 

• Moreover, the equation ()84p defines a formal Fourier series of a generating function which 
can be written as T" = Tq + T" where Tq has been defined in (|69p and Ti satisfies 

\KO^T^l^.,<Kl^G^\ 

for any < m + n < 1 . 

Proof. The first part follows from rewriting Lemma 15.191 in this setting and taking into account 
that ||F"||o,K,o- ^ KfiG^'^. The second, from the uniqueness statements of the first part and 
Theorem 15.161 

Finally, the third statement can be easily deduced from equation (j84p and taking into ac- 
count that the formal generating function was already defined in D^^ n -D^^^jj thanks to 
Theorem [5T01 ' ' ' □ 

With this proposition we have finished the extension procedure and we have both invariant 
manifolds parameterized as formal Fourier series in the domain 1)^2,(52 • summarize this fact 
in the following theorem. 

Theorem 5.22. Let k,2 and 82 the constants given by Proposition \5.21\ Then, there exist formal 
Fourier series T"''^ with Fourier coefficients defined in -0^2,52 which are solutions of equation ([24 
and satisfy 

\\^':^lT^'%,„,<h^^,G^\ 

with {) < m + n < 1 and > a constant independent of ji and Gq. 
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Proof. It is enough to join the results of Theorem 15.101 and Proposition l5.211 Indeed, in Theorem 
15.101 we have obtained solutions as formal Fourier series of equation (|24p of the form T"'* = 
Qq''^ + Q^'^ in the domains D^f^^ Using the fact that -Dk2,(52 <^ -^pi kq So estimates 
obtained in Lemma 15.131 and Theorem 15.101 we obtain the desired estimates for Tf . For the 
unstable manifold, we have to recall that -Dft2,52 ko 5o ^ ^k2,<52- Then, we need also to take 

into account the results obtained in Proposition 15.211 and recall that since the points -0^2,(52 
far from the singularities v = ±i/3, the norms || • ||i_o-2 and || • ||o,(T2 equivalent and satisfy 

II • ||o,o-2 ^ -^11 ■ lli,o-2- n 



6 The difference between the manifolds 

Once we have obtained the parameterization of the invariant manifolds (as formal Fourier series) 
up to points 0{Gq^) close to the singularities v = =ti/3, the next step is to study their difference. 
It suffices to study the difference between the Fourier series of the generating functions T"'* 
defined in ([23|) . To this end, we define 

A{v,0 = T'{v,0-T''{v,0- (85) 

Recall that T"'^(t;,,^) = To{v,^,) + T^'''{v,S,) where Tq is defined in ([69]) and T^'^ are the Fourier 
series obtained in Theorem 15.221 Therefore, the Fourier coefficients of A are defined in ,52. 
Recall that A is not a function but a formal Fourier series that does not need to be convergent. 

Subtracting equation for Tf{v,^) and T^{v,^) considered as equations of Fourier series, 
we can easily see that 

A G Ker£ 

where C is the differential operator 

£ = {1 + Aiv, 0) -Gl{l + B{v, 0) (86) 

with 

Biv,0 =^ ({d.Tl + d^T^) - ^{d^Tl + a^Ti")) (87) 



g; 



-3 



2r^ 



{d^Tf+d^Tr) 



where T^'^ are the functions given by Theorem 15.221 

Recall that the equation CA = is stated as an equation for not necessarily convergent 
Fourier series and not as an equation for functions. The derivatives for Fourier series have been 
defined in the natural way in ()63p . 

To study the function A we proceed as in |Sau01j (see also [GOSlOj ). That is, we consider 
a near the identity change of coordinates which conjugates the operator £ in (|86p with the 
operator C in ()34p . Recall that, as we observed in Section [531 the composition of formal Fourier 
series with a near the identity formal transformation is well defined under suitable hypotheses 
(see Lemma 16.21 below) . 

The structure of this section goes as follows. First, in Section [6. 11 we introduce a functional 
setting to study the difference between the generating functions. This functional setting is 
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essentially the same as the one considered in Section 15.4.11 but referred to Fourier series with 
coefficients defined in D^^^- In Section [6.21 we straighten the operator C. Finally, in Section! 
we use this result to prove Theorem 13.21 



6.1 Weighted Fourier norms and Banach spaces 

We devote this section to define Banach spaces for Fourier series with coefficients defined in 
Dk,(5- First, we define the Banach spaces for the Fourier coefficients as 

,i/^,K,5 — "[/i : Df^^^ y C : analytic, 

where 

\\h\\u.,u+= sup \{v-i/3r+{v + i/3r-h{v)\. 



Note that these definitions are the same as (|60|) and (j59|) but for functions defined in D^^a instead 
of ^ ^ . Now we define the Banach space for Fourier series 



+t/2,u-t/2,K,&^ \\ll'\\u,a < OO 

where 

The Banach space Qj/,k,(5,(t satisfies the algebra properties stated in Lemma 15.91 Therefore, from 
now on in this section, we will refer to this lemma understanding the properties stated in it 
as properties referred to elements of Q,v^n,&,a instead of elements of 3^i/,k,<5,ct- Moreover, in the 
present section we will need to take derivatives of and compose Fourier series. To this end we 
state the following two technical lemmas, which are equivalent Lemmas 15.141 and 15.151 

Lemma 6.1. Fix constants a' < a, k' > k and 5' > 6 and take h G Qu,K,s,a- Its derivatives, as 
defined in (|63|) . satisfy 



d'^h G Quy,S',a' and 



dfh G Qjy,K' 5' a' and 



1 



Lemma 6.2. We define the formal composition of formal Fourier series 

OO ^ 

h{v + g{v,0,0 = E —,d';:'Hv,Og^{v,0- 
^-^ ml 

m=0 

Fix constants a' < a, k' > k and 6' > 6. Let k' — k > rj > 0. Then, 
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If h £ Qu,n,s,cT, g G Qq,k',5',u' and ||g||o,cr' < tiGq^ we have that X{v,i) = h{v + g{v,C),£,) 
satisfies X £ yu,K',5',a' and 



X < - 1 



^/x./ „ x-i 



Moreover, if \\gi\\o,a, ||52||o,a < vGq^, then Y{v,C) = h{v + g2{v , ^) , - h{v + gi{v,C),0 
satisfies 

\\Y\\.y < (-Y (l - ' \\h\\u,.\\g2 - gi\\o,a. 



• If dyh G Q„,^,,5,cT, 51,52 G Qo,K.' ,5' ,ct' and 1151110,(7', 1152110,(7' < iiGq^ we have that Y G 
yu,K',S',a' and 

||i^||i/,(7' < (— j . rj \\dvHu,a\\g2 - gi\\o,a' ■ 

6.2 Straightening the operator C 

Once we have defined the Banach space Qu,K,s,a we can show that the operator C can be straight- 
ened, as stated in the foUowing theorem. 

Theorem 6.3. Let a2, k,2 and 82 he the constants given by Theorem \5.22l Let 1T3 < a2, H3 > K2 
and 5^ > 82 be fixed. Then, for Gq big enough, there exists a (not necessarily convergent) Fourier 
series C G Qo,k3,<53,<73 satisfying 

||C||o,(73 < bGfiG^\ 
with bQ > a constant independent of fi and Gq, such that 

A(u;,0 = A(u;+C(zz;,0,6, 



where A is the function defined in ()85p , is well defined and satisfies that A G Ker£ and C is the 
operator defined in ([31]) . 

Moreover, one can choose C to satisfy 

Ci-w,-0 = -C{w,0, 

and therefore, by ()27p . one has that 

A{-w,-0 = A{w,0- 

To prove this theorem we need a change of variables v = w + C{w, ^) such that A G Keri2 if 
and only if A G Keri2. This fact is equivalent to look for a Fourier series C which is solution of 
the equation 

£(C) = /C(C) (89) 

where 

r(hV A{v,0-B{v,0 



l + B{v,0 

and A and B are given in ([HI 
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We devote the rest of this section to obtain a solution of this equation. Let us point out that 
the fraction in the right hand side, since it involves formal Fourier series, is understood as 

The solution of equation (|89|) is found through a fixed point argument. As we will see in 
Lemma 16.61 the functions A and B inherit the symmetry properties of the manifolds T" and 
in (|27|) . To obtain a symmetric solution C of the equation (|89|) . we need a suitable left inverse 
of the operator C in the domain Dk,<5) which is a slight modification of the operator Q in ()66p 
acting over the Fourier coefficients and is defined as 

m{v,0=Y.^ihf^{v)e''^, (90) 

where its Fourier coefficients are given by 

g{hf\v) = r e^^^o{«-t)/iM(t) dt for £ < 



-I fV2 1 fV2 

G{hf\v)= h^'^\t)dt-- h^^]{t)dt-- h^°\t)dt 



2 J,* ^ ' 2 

V 



g{hf\v)= / e'^^o(^-'^h^i]{t)dt for£>0. 



V2 



Here V2 = ^(1/3 — kGq ^) is the top vertex of the domain Df^ s, V2 is its conjugate, which 
corresponds to the bottom vertex of the domain D^^s and v* is the left endpoint of D^^^s H M. 

Lemma 6.4. The operator Q in (|9U|) satisfies that if h & Qu,K,S,a for some v € (0, 1), then 
Qih) G Qo,K,5,(7, and 



Gih) 



0,(T 



< K\\h\L^. 



Moreover, 

• If h is a real-analytic Fourier series, that is 



then so is Q{h). 
• If h satisfies — ^) = h{v,^), one has that 

g{h){-v,-0 = -§{h){v,c). 

Proof. The first part is proven as Lemma 8.3 of jGOSlOj . The real-analyticity property is 
straightforward. To prove the symmetry property the first observation is that, if h is even and 
taking into account that V2 = —V2, one has that for ^ 7^ 0, 

g{h)^-'\-v) = -g{hf\v). 

Therefore, one just needs to see that g{h)^^\v) is an odd function. To this end, it is enough to 
check that 

f{v) = g{hf^{-v) + g{hf^{v) 

satisfies, wherever it makes sense, f'{v) = and f{v2) = 0. □ 
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We look for a fixed point of tlie operator 

K. = Q oK. 

in tlie space Qo,k,<5,o-- Tlieorem 16.31 is a straightforward consequence of next proposition. 

Proposition 6.5. Fix as < a"2, K3 > K2 o,nd 5^> 82- Then, there exists a constant h%> Q such 
that for Go big enough the operator K, is well defined from 3(1^^x0^'^) C Qo.ks.^s.o-s to itself and 
it is contractive. Therefore, it has a unique fixed point C in this ball, which is a real-analytic 
Fourier series and satisfies the symmetry condition C{—v, — ^) = — C(f,^). 

To prove the proposition, we first state the following technical lemma, whose proof follows 
easily from Theorem 15.221 and the properties of the homoclinic parameterization given in Corol- 
lary US] and (fT6]l . 

Lemma 6.6. The functions A and B defined in ([87j) satisfy that A G Qi/2,k2, 52,0-2 '^'^^ ^ ^ 
Qi/2,K2,<52,o-2 "^^^ symmetry properties 

A{-v,-0 = A{v,(,) and B{-v,-(,) = B{v,(,). 

Moreover, 

PlliAa2 < KfiG^' 

1 1 1 1 1/2,(72 ^ K^iGq^. 

Proof of Proposition 1 6. 5t First we see that if we fix any K > 0, then for Go big enough, 
the operator /C is well defined from B{K ^Gq^) C Qo,ki,<5,o-i to Qo,k3, 53,0-3- Indeed, take h G 
K3,53,o-3- Then, applying Lemmas [6^2] and [6^ with rj — K^Gq ^ and taking Gq 
large enough, one can easily see that K.{h) G Qi/2,k3,53,o-3- Then, applying Lemma [5. Ill we have 
that K,{h) = g o }C{h) G Qo,ki,5,oi- 

Once we know that the operator is well defined, we prove that it is contractive in a certain 
ball. As a first step we look for bounds for /C(0). By Lemma 16.61 one can easily see that 

Il^(0)||l/2,<X3 <^/^Go^ 

Then, applying Lemma [631 we have that there exists a constant b^ > such that 

11^(0) ||o,.3 < l/^Go^ 

Finally, we take /ii,/i2 G B{bQfiGQ^) C Qo,ft3,53,o-3 and we prove that K. is contractive. To this 
end, we first bound /C(/i2) — fC{hi). Using Lemma [6. 21 and the estimates obtained in Lemma [6.6^ 
one can easily see that 

||/C(/i2) - /C(/li)||i/2,.3 < K^iG^^\\h2 - /ll||o,<X3- 

Therefore, applying Lemma 16.41 we obtain that 

||/C(/l2) - /C(/li)||o,a3 < i^AiCg ^||/l2 - ^l||o,<73- 

Thus, this gives the existence of a unique fixed point C of the operator /C. Furthermore, by 
Lemmas l6.4l and [6r6| one has that IC(h)(—v. — f;) = —lC{h){v,^) if /i satisfies —^) = —h{v,^) 
and therefore the fixed point C satisfies the required symmetry property. Reasoning analogously 
one also obtains real-analyticity. This completes the proof of the proposition. □ 
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6.3 Estimates for the difference between invariant manifolds 



We devote this section to complete the proof of Theorem [221 that is, to prove that the generahzed 
Poincare function ()28p gives the first order of the difference between the manifolds. As a first 
step, we show that not necessarily convergent Fourier series ^' G Qo,K,s,a such that ^ G KerC 
define functions for real values of the variables which have exponentially small bounds in Gq. 

Lemma 6.7. Fix k > 0, 5 > and a > 0. Let us consider a formal Fourier series ^ £ Qo,k,s.ct 
such that ^ S Ker£. Then, the Fourier series ^(i),^) 



is of the form 



for certain constants G C and 

defines a function for v G -0^.(5 H M and G T, which satisfies that 



< sup 

< sup 



< 



sup 



i^l^lG^e 3 



Proof The fact that ^ G Ker£ implies that each Fourier coefficient ^'[^^ satisfies 

d^^^^^v) - UGl^^^v) = 
and therefore, there exist constants A^ such that 

Evaluating this equality at the top vertex V2 = ^(1/3 — kGq^) of -Dk,5 for £ < and at the 
bottom vertex V2 = — i(l/3 — kGq^) for £ > 0, we obtain that 



This implies that for v G -D^ 5 H ^ 



aM 



e ' ' 3 



< sup 



Moreover, since ^' G Qo,k,<5,(ti we have that 



From this bound one has that, for v G H M and G T, the Fourier series of ^ is convergent 
and therefore it defines an analytic function. Finally, taking derivatives one can easily prove the 
bounds for 9^^-^ and d^¥^l □ 
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We use Lemma 16.71 to prove that ()28p gives the first order of the difference between the 
manifolds. First let us observe that, since by Theorem \<j.'6\ the Fourier series A in (|88p is 
symmetric, real-analytic and satisfies A G KerC, it is of the form 

A{w, = r[°i + 2 ^ rW cos e {gIw + ^) . (91) 

for certain coefficients G M. 

A direct application of Lemma [6.8| provides exponentially small bounds for A. Nevertheless, 
our goal is not give bounds but to prove that the function L in ()28p is the main term in A. 
Thus, we define the formal Fourier series 

Note that L can be expressed as a formal Fourier series with coefficients in the domain D^ s ^ 

L = {Qt>-Q^)-{Ll-^), (92) 

where Qq and are the Fourier series defined in (|39p and (j52p . respectively, and Qq and Lf 
are the analogous Fourier series referred to the stable manifold. 

Lemma 6.8. Consider the constants K3 and 63 defined in Theorem 1 6. 3i Then, for {w,^) G 
K3,53 nM) X T, the Fourier series £ defines an analytic function, which satisfies 

\£{w, 0-E\< Kfi^ (1 - 2^) Go ^e-^ + KGq ^^V^e"^, 



where E &W is a constant, and for < m + n < 2, 



3 



Proof. To proof the lemma it is enough to point out that £ G Ker£ and then apply Lemma 16. 71 
We use the expression of L given in (|92p . By Theorem 16. 3t we have that A G Ker£. Moreover, 
C{Qq) = Uq for both * = u and * = s and therefore C{Qq — Qq) = 0. Reasoning analogously, 
C{Lf — Lf) = 0. Then, to apply Lemma 16.71 it only suffices to bound ||<?||o,o-3 in the domain 
I?K3,<53 n ^p2,K3,53- To this end, we split £ as 

£ = £(- £^ + £2 

where 

£1 = Tl -Ql-L\ 
£2 = A- A. 

We start by bounding ||iS"||i^o-3 • ^y the definition of in (P0|) . we have that = — L". 
Then, applying formula ([5^. we obtain 

iifriii,.3 <^yGo^ 

Then, applying Lemma 15.91 we obtain that 
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The bound for is analogous. To bound £2, recall that by the definition of A in 
be written as 

S2{w, = ^{w + C{w, 0,0 - ^{w, 0- 



it can 



Prom Lemma [5TT3I Theorem 15.101 and Lemma [521 we know that d^A G Q; 



3/2, ft3 ,53,0-3 



and satisfies 



||9„A||3/2,o-3 < K^lGq^. Therefore, applying Lemma [62] and using the estimates for C given in 
Theorem 16.31 we obtain that 



Il'^2||3/2,a3 < 

Applying Lemma 15.91 we obtain that 



13/2,(72 



|C||o,.3 < Ki?G^\ 



and thus we can conclude that 
which implies 



\£\\o,a,<K^?G^'"\ 



7 I 3 
2 



Recall that, when /u = 1/2, the Hamiltonian ()17p is vr-periodic. It can be easily seen that this 
fact implies that £ is vr-periodic in ^ and therefore, when fi = 1/2, it satisfies that £^^^ = for 
any odd £. On the other hand, the function £ depends analytically on ^ and therefore one can 
apply the Schwarz Lemma to see that for any odd i, 



l+fl 



£^^^{v) < fi^{l-2fi){KGo 
Applying now Lemma \67l\ we obtain that for w G Df^^^s-s H M, 



£^'Kv) <^'(l-2^f)(KGo)-i+5l^le- 



< f^^KGo) 



7 , 3|/,| _KKlO 

'2 + 2l*le 3 



for odd 



for even \i\ 



and analogous bounds for the derivatives of the Fourier coefficients. Defining E = and 
summing up the odd and even Fourier coefficients of £ one obtains the bounds stated in 
Lemma 16. 8[ □ 



Now it only remains to go back to the original variables {v,0- This is summarized in the 
next lemma, from which the proof of Theorem 13.21 follows. Recall that now we are taking real 
values of the variables and therefore all the objects we are dealing with are functions and not 
only formal Fourier series. 



Lemma 6.9. Consider the functions A defined in (|85|) and L defined in (|28|) and the constants 
K3 and 63 given by Theorem 1 6. 31 Fix K4 > K3 and 64 > 5^. Then, for {v,0 £ (-Dk4,(54 H M) x T, 



l/2„2„-^ 



A{v,0 - L{v,0 - E <Kfi^{l-2fj,)GQ^e-~+KG(^ ' fi'e 
for certain £' G R, and for < m + n < 2, 

d'^dlA{v,0 - a™6>^"L(t;,0 < Kf? (1 - 2^l) Go 2+3^6-^ + J^Go'/'+^'^/i^e-^, 
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Proof. It is enough to consider the inverse of the change of coordinates (f , ^) = {w + C{w,(,),S,) 
obtained in Theorem 16.31 Recall that now we are interested in real values of the variables. 
Nevertheless, since we need to apply Cauchy estimates we consider a small complex neighborhood 
of (D^a.^s nM) X T. Then, shrinking slightly the domain, one can easily obtain the inverse change 
by means of a fixed point argument. It is of the form (w,^) = {v + T{v,(,),^) and T satisfies 
1 1 I loo ^ K^iGq^. Applying Cauchy estimates we also know that for real values of the variables 
ll^ioTf^lloo < K^Gq^ and ||c?^T||oo < K^Gq"^. Applying the change of coordinates to the bounds 
obtained in Lemma 16.81 and using the bounds obtained for T and its derivatives, Lemma 16.91 
easily follows. □ 

Proof of Theorem \2.S[. Due to the form of A in (|9ip , one has that 

A{w,i)=T[Glw + ^) 

Moreover, it is clear that T'['Kk) = V"'{Tik) = for any A; G Z. On the other hand, by Lemma [6.9l 
and Proposition 13.1 1 on has that 

T"{z) =2^fi{l - fi)G~^^^e-^ 
+ 0(MGfe-«g). 

Therefore, for any odd A; G Z, the equation T"{kTT) = defines a curve t] in the parameter plane 
(//, Go), which is of the form 

A^ = f^*iGo) = \- 16V2Gge-# {l + O (g^^)) . 

Moreover, for (^, Go) € r/, one has that r(™)(A;7r) ^ 0. This implies that all partial derivatives 
up to order three of the function A are zero and the fourth are not whenever GqW + ^ = /ctt for 
odd k. 

Now, for Go big enough, we can choose k such that, for any ^ G T, one has that = 
Gq ^(/cTT — .^) G -Ck3,<53 n M . Therefore, taking into account the relation between A and A in 
([88]) . at the point {vk,C) = {wk +C{wk,£,),C), all the partial derivatives up to order three of the 
function A vanish and the fourth one do not. 

To finish the proof of Theorem 12.31 it only remains to go back to the parameteritzations (|19p 
of the stable and unstable invariant curves 7"'* of the Poincare map (jl8p . They are defined by 
formula ()29p . Then, at the point v = Vk, one can easily see that 

{Yl{v,;fi*{Go),Go)-Y^^ivk;fi*{Go),Go)) = 
92 {Yl{vk;fi*{Go),Go) - Yl{v,; fi*{Go),Go)) = 
{Yl{v,;fi*{Go),Go)-Ylivk;fi*{Go),Go)) / 0. 

Therefore, the invariant curves 7"'** have a cubic homoclinic tangency. 

Let us emphasize that, for (/i,Go) G rj and any even k, T"{k'Tr) 7^ 0. This implies that, at 
the corresponding homoclinic points, the invariant curves 7"'^ intersect transversally. □ 



1-2^ 



cos , 



1 + O G, 



,-1/2 

'0 



+ 8G^ 



3 



cos 2z 1 + O G, 



-1/2 
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A Computation of the function L: proof of Proposition 13.1 



We devote this appendix to give the properties of the function L in stated in Proposition l3.1l 
We first state an auxihary lemma, which gives some properties of the potential V in (jlip . 

Lemma A.l. The Fourier coefficients of the function 
where V is defined in (jlip and Fh in ()10p , are 

w/iere Cj = ^ "^^ ^ , So{0) = 1 and 6o{£) = Ofori^O. 

Proof. We use the identity 

(1 + Acos = «-i/2(i + /3e^e)-i/2(i + /3e-''^)-i/2 



a 



j>0 A;>0 



j>max{0,-^} 



where a = A/ {2(5) and /? = (1 - Vl -^2)/^, and 



V^G4f2_2^G^fcose + /i2 Ggr V Ggr y V ^gr 



VGJF^ + 2(1 - f,)GJFcos9 + (1 - m)2 V G^r y V G^r 

to obtain 



V^G4f2_2^Ggfcose + /z2 v^G^fS + 2(1 - /i)G2~pQg ^ + (l _ ^)2 G2r(?;) 

□ 
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Now we prove Proposition 13.11 
Proof Proposition \3. 1[ First, we observe that the function L in ()28p can be written as 

/+00 
y(fh(t), e + Glv - Git + ah{t);n, Go)dt. 
-oo 

Now, using that 

V{¥Ut),^ + Glv - Git + 5h(0; /i, Go) = C/M(t)e^^"'^(*)e-^^^o*e^^«+Gg-) 
and therefore 
with 

r+oo 



to compute we use the expansions in Lemma lA.ll obtaining, when i ^ 
with 

To compute we use the change of variables v = ^ (5''"'^ + ''") given in Lemma [4.11 obtaining 

-^-,-1— i™dr := (-l)^22^+^/(-Ai,J +^), 



where we have introduced the notation 



I(i, m,n) = / -75 — ; r^:7^dT 



to write the Melnikov potential Fourier coefficient as 



j>0 ^ ^ ^ - ^ ^0 
The first observation is that 



/(— n, m) = m, n) = m, n). 
Therefore all the Fourier coefficients are real and we just need to compute them for £ > 0: 
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and then 

L(^, z;) = 2 ^ lM cos + Glv) . 

im 

To compute the integrals m, n) for ^ > 0, one uses the method in |Erd56] (see also 
|LS80al [MP94j ) changing the path of integration to a suitable complex path Re (r + ^) = up 
to a neighborhood of the singularity t = i. 

Using that t + ^ = + 0{{t - i)'^), to bound the integrals (see |DKdm,S12j ^ it is enough 

to reach a neighborhood of the singularity r = i of order 0{Gq ), obtaining that there exists 
a constant K > such that, for any ^ G Z and m,n > 1: 



\Iii,m,n)\ < KGo^-'^^^e-^ 

and therefore 



g: 



3 



< {KGo 



3 



To obtain the dominant terms of the function L, which correspond to £ = 1,2, we can use the 
results in |DKdlR,S12j to obtain 



1 In _l ''o 



1(1,2,1) = l^/l-Gle-^ (l + 0(Go"'/' 



6 V 2 



/(2, 2, 0) = 2V^G^/'e-2% (l + O (g" 



Thus, 



L^l = 2/z(l - /x)V^G|e-2% (l + O (Gq 



3/2 



□ 
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